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Abstract: 

The elliptic Einstein-DeTurck equation may be used to numerically find Einstein metrics 
on Riemannian manifolds. Static Lorentzian Einstein metrics are considered by analyti- 
cally continuing to Euclidean time. Ricci-DeTurck flow is a constructive algorithm to solve 
this equation, and is simple to implement when the solution is a stable fixed point, the 
only complication being that Ricci solitons may exist which are not Einstein. Here we 
extend previous work to consider the Einstein-DeTurck equation for Riemannian mani- 
folds with boundaries, and those that continue to static Lorentzian spacetimes which are 
asymptotically flat, Kaluza-Klein, locally AdS or have extremal horizons. Using a maxi- 
mum principle we prove that Ricci solitons do not exist in these cases and so any solution 
is Einstein. We also argue that Ricci-DeTurck flow preserves these classes of manifolds. As 
an example we simulate Ricci-DeTurck flow for a manifold with asymptotics relevant for 
AdS^/CFT^. Our maximum principle dictates there are no soliton solutions, and we give 
strong numerical evidence that there exists a stable fixed point of the flow which continues 
to a smooth static Lorentzian Einstein metric. Our asymptotics are such that this describes 
the classical gravity dual relevant for the CFT on a Schwarzschild background in either the 
Unruh or Boulware vacua. It determines the leading 0{N^) part of the CFT stress tensor, 
which interestingly is regular on both the future and past Schwarzschild horizons. 
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1. Introduction 



Exact solutions to Einstein's equation are notoriously difHcult to find and traditionally 
most progress has been made under extra symmetry assumptions. However, in recent 
times it has become apparent that there are many interesting solutions of the Einstein 
equations with less symmetry than the usual examples. -"^ Wc arc likely never to know these 
analytically, and so if we wish to have a detailed understanding of such - believed to exist 
- solutions, one must resort to numerical techniques. In modern times there are three key 
areas where such questions arise: i) string theory phenomenology, where one is interested 
in the Riemannian geometry of the extra dimensions, ii) black holes in higher dimensions, 
where there is a very complicated phase structure, even in the static case if the extra 
dimensions are compact, and iii) AdS-CFT, where exotic gravitational solutions of many 
types have gained interesting physical significance due to their relation to certain strongly 
coupled CFTs. 

A general approach for constructing static vacuum solutions has been outlined in [1]. 
The static spacetime is analytically continued to a Riemannian geometry and then the 
elliptic Einstein-DeTurck equation is solved as a boundary value problem. This approach 
is particularly attractive because if the spacetime contains a non-extremal black hole, after 
analytic continuation with suitable Euclidean time period, there is no boundary associated 
to the horizon location and the Riemannian manifold is smooth there. Two numerical 
algorithms have been proposed to solve the Einstein-DeTurck equations: simulation of the 
Ricci-DeTurck flow, and the Newton method. The Ricci-DeTurck flow is the most elegant 
since it is geometrical, and turns out to be simple to simulate numerically. It is this method 
that we focus on in this paper. A number of previous works have computed Ricci flows 
numerically in a variety of related contexts [1-7]. 

One aim of the current paper is to extend the discussion of the Einstein-DeTurck equa- 
tion in this numerical context to allow for new types of boundary or asymptotic conditions, 
beyond the simple one considered in [1]. We will restrict to vacuum solutions with zero 
or negative cosmological constant. We show how to impose boundary conditions for gen- 
eral asymptotically Euclidean, Kaluza-Klein and locally hyperbolic Riemannian manifolds 
which continue to give Lorentzian geometries which are asymptotically flat, Kaluza-Klein or 
locally AdS. We also show how to deflne boundary conditions for static Riemannian mani- 
folds^ whose Lorentzian analytic continuations contain extremal horizons (this corresponds 
to an asymptotic region in Euclidean signature). We also consider the Einstein-DeTurck 
equation on general manifolds with boundaries. Wc discuss the boundary conditions con- 
sidered by Anderson [8], where the conformal class of the induced metric and trace of 
extrinsic curvature are fixed. We also examine the case where the extrinsic curvature is 
proportional to the induced metric. 

Furthermore, since we wish to use Ricci-DeTurck flow to solve the Einstein-DeTurck 

^The simplest example of this is five dimensional KK static vacuum black holes, where on top of the 
uniform black string, one can have non-uniform black strings and localised black holes. 

^We will refer to a 'static' Riemannian geometry as a Riemannian manifold possessing a 1/(1) or R 
isometry generated by a hypersurface orthogonal Killing field. 
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equation, we have considered whether our boundary conditions are preserved under such 
flows. For the asymptoticaUy Euchdean case this has been considered previously in the 
Uterature [9] where it was shown that Ricci flow preserves this class on manifolds. It seems 
reasonable to expect that this can be generalised to the Kaluza-Klein case, although we have 
not considered this. The asymptotically locally hyperbolic class has recently been proven 
to be preserved by Ricci flow under a smoothness assumption [10]. For static manifolds 
with asymptotic regions that analytically continue to Lorentzian extremal horizons, we 
have proved here that they are indeed preserved by Ricci-DeTurck flow; the key idea here 
is that there exists a well defined general notion of a near-horizon geometry which itself 
solves the Einstein equations [11]. 

There arc two subtleties to using Ricci-DeTurck fiow to solve the Einstein-DeTurck 
equation. Firstly a static solution in vacuum gravity may not be a stable fixed point of 
the flow.^ This is often the case for black hole solutions due to the existence of Euclidean 
negative modes [12]. In [1] it was proposed that by appropriate tuning of a family of initial 
data one may still find solutions, although in practice the Newton method usually works 
best in these situations. Therefore Ricci-DeTurck flow is best applied to cases where the 
fixed point of interest is believed to be stable. 

A second subtlety is that there exist solutions of the Einstein-DcTurck equation that 
do not correspond to solutions to the Einstein equations, but instead arc Ricci solitons. 
To address this, we show that there is a simple maximum principle which impacts on the 
existence of Ricci solitons. Applying this maximum principle to the cases with boundary 
and asymptotic conditions considered above, in fact rules out the existence of Ricci solitons. 

In order to illustrate these general ideas, we consider an interesting example of a 
Ricci-DeTurck flow. Specifically it is a flow on a five dimensional static, axisymmetric 
Riemannian manifold with two asymptotic regions that is smooth in the interior. One 
asymptotic region is locally hyperbolic with Schwarzschild boundary metric, and the other 
has a Lorentzian analytic continuation which gives an extremal horizon whose near-horizon 
geometry is that of standard Poincare AdS. By our general results we know how to define 
boundary conditions such that the flow preserves such a class of Riemannian manifolds. 
Furthermore, from our non-existence result concerning Ricci solitons, we know that if the 
flow converges to a fixed point, it must be an Einstein metric. We provide strong numerical 
evidence that the flow does converge, and that such a fixed point with smooth interior does 
exist. Furthermore, by the above comments we expect this fixed point is stable, since no 
fine tuning was required to flow to it. 

The Lorentzian continuation of this new Einstein metric has a very interesting AdS^/CFT4^ 
interpretation. Namely, as the regular static classical bulk geometry dual to the Lorentzian 
signature CFT4 on a Schwarzschild black hole background in a vacuum state. Indeed, the 
bulk geometry has an asymptotic AdS boundary in the 'UV end of the geometry, with 
conformal boundary metric being Schwarzschild. Furthermore, since in the bulk the geom- 
etry approaches that of the Poincare- AdS horizon, it is conformal in the TR' and hence is 
dual to the CFT in a (non-thermal) vacuum state. We argue that this classical solution 

^Recall that the stabihty (under the flow) of a given fixed point is determined by the non-existence of 
negative eigenvalues in the spectrum of the Euclidean Lichnerowicz operator for that fixed point. 
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describes the dominant bulk saddle point in AdS^/CFT^ suitable for describing the Unruh 
and Boulware vacua in the CFT. Since the solution has an extremal horizon in the IR, we 
do not expect it gives the dominant saddle point for the CFT in the Hartle-Hawking vac- 
uum, which presumably would be of the form of a 'funnel' or 'droplet' discussed in [13, 14] 
with a non-extremal horizon in the IR with the same temperature as the boundary black 
hole. 

A classical bulk solution allows one to compute the leading 0{N^) part of the dual 
CFT stress tensor, and we do so from our numerical solution. We believe this is the 
first calculation of the leading order contribution to a stress tensor at large Nc for a four 
dimensional strongly coupled field theory in an asymptotically flat black hole background. 
Since our bulk solution is smooth in the interior, we find the remarkable result that this 
leading contribution to the stress tensor is regular on both the future and past horizons even 
though the solution describes the dominant classical saddle point for the dual to the Unruh 
and Boulware vacua. Since our calculation gives only the leading O(iV^) contribution at 
large Nc and is at strong 't Hooft coupling, the familiar free field theory properties of the 
Unruh and Boulware vacua are not apparent (i.e. the former being regular only on the 
future horizon and the latter regular on neither future nor past horizon). 

Finally, we emphasize that this example represents an optimal situation for the nu- 
merical construction of a new solution which is of interest physically, extremely simple to 
find just using Ricci-DeTurck flow (as it is a stable fixed point), and where one can rule 
out the potential complication of existence of Ricci soliton solutions. 

This paper is divided into two main parts. In the first part we extend the analysis of [1] 
to understand how to impose the various boundary conditions of interest for the Einstein- 
DeTurck equation, and for the Ricci-DeTurck flow. This begins in §2.1 where we point 
out a simple maximum principle for the Ricci soliton equation. Then in §2.2 we consider 
various new boundary and asymptotic behaviours for the Einstein-DeTurck equation, and 
show that the maximum principle applied in such cases rules out the existence of soliton 
solutions of the Einstein-DeTurck equation. Then in §3 we consider Ricci-DeTurck flow as 
a method for solving the Einstein-DeTurck equation, and discuss whether our boundary 
conditions are preserved by such flows, proving this is the case for static manifolds whose 
Lorentzian section contain extremal horizons. In the second part, §4, we illustrate these 
general ideas by giving an example of a Ricci-DeTurck flow, and provide strong numerical 
evidence that it converges to an Einstein metric with both asymptotically extremal and 
locally hyperbolic regions. The Lorentzian continuation of this is a smooth Einstein metric 
which is asymptotically locally AdS with Schwarzschild boundary metric and also contains 
an extremal Poincare-AdS horizon. We conclude the paper in §5 with a brief summary 
and a discussion mainly devoted to considering the AdS/CFT interpretation of our new 
solution. Some of the technical details are relegated to several appendices. 

2. Einstein-DeTurck solitons 

2.1 A maximum principle for Ricci solitons 

Let {Ai,g) be a smooth (connected) m-dimensional Riemannian manifold which obeys the 
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Ricci soliton equation 



(2.1) 



for some 1-form ^, where V is the Levi-Civita connection associated to the metric g. We 
will consider a cosmological term which is either vanishing or negative (hyperbolic) , so that 
A < 0. We term a solution (g,^) of this equation with non-vanishing ^ to be a non-trivial 
soliton. 

The contracted Bianchi identity implies that any solution to the soliton equation must 
satisfy the vector equation, V^^^ + R^£,i, = 0. Combining this equation with the original 
Ricci soliton equation (2.1), and defining the function cj) = ^°^^a ^ 0, one can establish the 
scalar equation 



where we have noted that both terms on the right hand side are non- negative. 

Clearly for a solution to be a non-trivial soliton it is necessary (although certainly not 
sufficient) for the associated equation 



for a function / on the fixed background (7W,g,^), to have a non- vanishing solution such 
that / > 0. Furthermore, the behaviour of / on any boundary dM, needs to be the same 
as that for cj). However, this equation admits a maximum principle^ which for non-constant 
/ states: i) / may attain its maximum only on the boundary dM, and ii) the outer normal 
gradient obeys 9„/ > at such a maximum on dM. Thus, in particular note that if / > 
and / = on dM, statement (i) implies that / = everywhere in M. 

On a compact manifold with no boundary, the maximum principle thus shows that / 
must in fact be a (possibly non-zero) constant, and therefore in the soliton case one learns 
that is a non-negative constant. Furthermore, if A < the r.h.s. of (2.2) in fact implies 
(j) = Q (and hence the soliton is trivial), whereas if A = one has V^^,/ = (and hence 
the soliton is Ricci fiat). In fact compact solitons with ,^ 7^ do not exist for A = in 
the absence of boundaries [17], although that isn't apparent from our simple maximum 
principle argument. 

Our interest here is in both manifolds with boundaries and non-compact manifolds. 
In particular, we are interested in general manifolds with a boundary, such that either 
(/> = on dAi, or the outward normal derivative (9.„(/) < on dAi. We are also interested 
in non-compact manifolds with a number of ends, such that (p ^ Q asymptotically in each 
end. In these cases, from above and as we argue below, the maximum principle rules out 
a non-zero solution to (2.3), and hence non-trivial soliton solutions. 

^Hopf's maximum principle is the relevant local result and can be found in standard texts [15]. For 
related global statements on Riemannian manifolds see [16]. 





+ ed^^f > 
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2.2 Einstein-DeTurck equation on manifolds with boundaries and asymptotic 
regions 

From now on we consider the DeTurck choice of vector field ^, which is constructed from 
the Levi-Civita connection T of g and a smooth reference connection F, 

e = g'-' {k, - f . (2.4) 

This was introduced by DeTurck in the context of an elliptic problem, and later used by 
him to argue Ricci flow is diffeomorphic to a parabolic flow (Ricci DeTurck flow) [18]. 
In this case we term the Ricci soliton equation (2.1) the Einstein-DeTurck equation, and 
unlike the usual Einstein equation, for Riemannian signature this is elliptic. Hence we 
expect to be able to solve it for suitable boundary conditions as a boundary value problem. 
We further make the choice that the reference connection is the Levi-Civita connection of 
a reference metric g on M. 

Before we give some examples of boundary and asymptotic conditions of interest, it is 
useful to introduce a general chart («;,x*) for the metric^ 

g = a^{dw + Pidx^f + -fijdx^dx^ (2.5) 

where the coordinate w controls the proximity to the boundary or asymptotic region. We 
may always choose coordinates such that the reference metric is in normal form, 

g = cP'dw^ + jijdx^dx^ . (2.6) 

Now working in the non-coordinate basis, 

e'^ = dw + Pidx'' , e* = dx* , eu, = du, , Si = di - Pidu, , (2.7) 

we can calculate, 

= - - - (2-8) 

^^ = 5^ A - + 1^ + I'^'i^ikl - ^^kl) (2.9) 

and 

</. = ieP = ^(^^)'+7'^e^e,. (2.10) 

In order to use our maximum principle we need (f) to vanish on the boundary, and for 
regular boundary metric this implies we will need and to vanish separately. 

2.2.1 Manifolds v^'itYi a boundary 

Consider the case where the manifold has a boundary dM. For the DeTurck choice of ^ 
the soliton equation is elliptic, and hence boundary conditions for all components of the 
metric g should be specifled on dM. 



^We emphasize that this is not the ADM decomposition of the metric, so in general 7ij is not the induced 
metric on the boundary and a, Pi are not the lapse and shift. 
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However, geometrically one considers only imposing boundary conditions on cither the 
induced metric or the extrinsic curvature. In D dimensions the full metric has D{D + 
l)/2 components, and the induced metric or extrinsic curvature each have only D{D — 
l)/2, the difference corresponding to the diffeomorphism freedom normal to the boundary. 
Thinking of an ADM decomposition with the normal direction to the boundary being 
'time', the lapse and shift are these non-geometric data which must be fixed in order to give 
enough boundary conditions for the elliptic Einstcin-DeTurck equation. These remaining 
D components are fixed by imposing boundary conditions on the ^ vector. We wish to find 
Einstein metrics and hence it is natural to require conditions on ^ at the boundary. For 
example, demanding ^ vanishes provides precisely the additional D boundary conditions 
needed to specify the lapse and shift. 

In the context of the Einstcin-DeTurck equation we will now consider two natural 
sets of boundary conditions that might loosely be termed 'modified Dirichlet' and 'mixed 
Dirichlet-Neumann' conditions. 

'Modified Dirichlet': The most obvious boundary condition would be to fix the induced 
metric, and also demand that C = 0- This locally provides the correct number of conditions 
for the metric components. However Anderson has shown [8, 19] the surprising result that 

fixing these 'Dirichlet' boundary conditions does not result in a regular elliptic problem for 
the Einstein equations modified by a Bianchi gauge fixing term, which is closely related to 
the DeTurck term we consider here when then metric is near to the reference metric.^ In 
fact more generally Anderson has argued that fixing the induced metric does not result in 
an elliptic problem for Einstein metrics, and an obstruction to ellipticity is given by the 
Hamiltonian constraint normal to the boundary. 

Anderson proposes that rather than fixing the entire induced metric, instead one can fix 
the conformal class of the induced metric, and the remaining freedom is fixed by specifying 
the trace of the extrinsic curvature. We term these 'modified Dirichlet' conditions. For 
the Bianchi gauge Einstein equations he has shown that together with requiring ^ = this 
gives a regular elliptic system. Following Anderson's analysis in [8] it is a simple exercise 
to confirm that the Einstein-DeTurck equation is also a regular elliptic system for these 
boundary data. 

We have therefore suitable elliptic boundary conditions for all metric components of 
g which precisely fix the conformal class of the induced boundary metric, the trace of the 
extrinsic curvature and, most importantly from our point of view, ensures that ^ = on 
the boundary. If the only boundary of M. is of this modified Dirichlet type, our maximum 
principle then implies there are no non-trivial solitons. If there are other types of bound- 
ary or asymptotic region, then these modified Dirichlet conditions ensure that if there is a 
non-trivial soliton, the maximum of (f) must reside on one of these other boundaries or in 
these asymptotic regions. 



We are grateful to Michael Anderson for pointing out an error in this subsection of a previous version 
of this work, and for explaining his work [8, 19] to us. 
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'Mixed Dirichlet-Neumann': We will now consider boundary conditions where we 
require the induced metric hij on the boundary to be proportional to the extrinsic curvature 

Kij + Xhij = 0. (2.11) 

Such boundary conditions arise in physical problems in the context of codimension one 
orbifold planes where the tension is related to A. An example of such boundaries arise in 
the Randall-Sundrum braneworld models [20,21] in the absense of matter on the branes. 
These orbifold branes may have both positive and negative A, corresponding to positive 
and negative tension. 

Let us first consider the case Kij = on the boundary. Whilst the manifold really ends 
at the boundary, in this case we may consider 'doubling' the manifold and smoothly gluing 
it across the boundary. This 'doubled' manifold no longer has a boundary but instead a 
reflection plane so that the metric has a discrete Z2 symmetry. This symmetry implies 
that the normal component of ^ vanishes, and the remaining tangential components of 
^ have vanishing normal gradient. Then, provided the reference metric shares the same 
discrete symmetry, the Einstein-DeTurck tensor on this doubled manifold will also haver 
this symmetry. With no boundary, the Einstein-DeTurck equation is then clearly elliptic. 
Furthermore applying our maximum principle to this 'doubled' manifold we deduce that 
(j) could not have a maximum on this reflection plane. Using the metric form (2.5) where 
■u; = is the reflection plane, we see that a and jij will be even functions of w, and /3i odd. 
Thus /3j vanishes at w = 0, and hence then the induced metric hij = 'jij on the boundary. 
Thus Kij = becomes d^'Jij = at w = 0. Likewise for the reference metric we require 
a and jij are even. We then explicitly see from (2.8) that = at w = and that is 
even in w. This example, Kij = at w = 0, is instructive as it shows that in this case, 
instead of imposing = on the boundary w = 0, it is more natural to fix ,^^„ = and 
Pi = at w = (which is consistent with the reflection symmetry of the doubled manifold 
introduced above). 

We now consider boundary conditions for our more general mixed Neumann-Dirichlet 
problem (2.11). We take our metric to be of the form (2.5) with the chart covering the 
manifold for w > and the boundary being u; = 0, so the outer unit normal 5„ = — ^9^- 
We begin by requiring that our reference metric (2.6) also obeys the same condition (2.11) 
as the metric, which implies Kij = —^d^jjij = —X^ij. We then impose the following 
boundary conditions on the metric (2.5) at u; = 0; 

= (aduia - ady^a) - A (a(m - 1) - aY^jij) = 
A = 

Kij + A hij = -^du,^ij + ^lij = (2-12) 

where m is the dimension of the manifold. We emphasize that since Pi = at w = 0, 
then the induced metric hij = ^ij and the extrinsic curvature Kij = —^d^jij at w = 0. 
Notice there are no tangential derivatives in these boundary conditions due to the condition 
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/3j = 0. In order to establish elliptic regularity of the system it is sufficient to consider the 
leading order derivative terms in the Einstein DeTurck and boundary operators acting on 
a, (3i and (the former of course acts diagonally). In other words, one must consider the 
short distance behaviour of the field equations and boundary conditions. Provided q 7^ 
(as required for a smooth boundary geometry), this yields, 

a^a~0, /3i~0, 9^7u-0 (2.13) 

where by ~ we mean the righthand side only contains terms without derivatives (which 
may be ignored at short distances). Then we see that linearised perturbations of this 
system satisfy Dirichlet conditions for the components of /3j and Neumann conditions for 
a and the components of ^ij . Therefore these boundary conditions give a regular elliptic 
boundary value problem. 

A straightforward calculation shows that at the boundary w = the tangential com- 
ponents of the vector ^ obey the simple condition, 

^n^^ = -2\ii (2.14) 

where to derive this expression one must use the boundary conditions above, together with 
the Einstein-DeTurck equations. Together with = these imply that, 

a„(/> = -2A0. (2.15) 

In the case that A = or A > 0, the zero or positive tension cases for Randall-Sundrum 
branes, then our maximum principle dictates that if a non-trivial soliton solution exists, 
the maximum of (p cannot reside on these boundaries, since the boundary conditions imply 
the outer normal gradient is zero or negative, whereas the maximum principle states at a 
maximum on a boundary the outer normal gradient must be strictly positive. If there are 
no other boundaries or asymptotic regions in the problem then this implies no non-trivial 
solitons can exist and hence any solution must be Einstein. If there are other boundaries 
or asymptotic regions this implies that if a non-trivial soliton exists, the maximum of (j) 
must reside on one of these other boundaries or in the asymptotic regions. 

We note that for the case of a negative tension orbifold plane, so that A < 0, then the 
boundary condition does not rule out the existence of a maximum of (j) on the boundary. 
Obviously (2.14) is compatible with finding Einstein metrics where = everywhere, but 
our maximum principle can not be used to rule out solitons in this case. In a numerical 
setting if a solution is found, one must explicitly check whether it is a soliton or not by 
computing (j). 

We note that Anderson has shown [8] that in the closely related Einstein-Bianchi 
equation the Neumann boundary condition together with ^ = is not regular. One can 
check this is the case also for the Einstein-DeTurck equation. We emphasize that for our 
boundary condition above we do not impose ^ = 0, but only the normal component S^yj = 0, 
and this is the reason we have regularity. The elegant behaviour (2.14) for the tangential 

^For a reference discussing regularity of PDEs see the textbook of Taylor [22] . 
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components of ^, which is consistent with finding Einstein solutions, requires our specific 
form of boundary condition, and presumably does not generalize to cases such as Kij = Sij 
for some tensor Sij on the boundary, unless Sij obeys special properties. Hence one cannot 
use our mixed Neumann-Dirichlet condition above to study the general Neumann problem 
with sources. 

2.2.2 Asymptotically Euclidean and Kaluza-Klein manifolds 

We now consider asymptotic behaviour of Riemannian manifolds that is appropriate to 
describe the Euclidean continuation of static asymptotically flat or Kaluza-Klein Lorentzian 
spacetimes, with or without a non-extremal horizon (so that Euclidean time may or may 
not be periodic). 

Consider the case where A = and (7W , (7) is a non-compact asymptotically Euclidean 
(AE) or asymptotically Kaluza-Klein (AKK) Riemannian manifold (where for the sake of 
generality we allow for a general compact Ricci flat internal manifold). By definition (see 
e.g. [23]), such a manifold M. has an end diffeomorphic to E = (W^"^ — Br) x Xj., where 
Bji is a ball of radius R and Xf; is a compact A;-manifold which admits a Ricci flat metric 
gx- For the AE case, k = and the end is simply diffeomorphic to E = (W" — Bji). 
Furthermore, if we let r = a/ Sijx'^x^ where are standard Cartesian coordinates on M"*~'^, 
the metric g in E, i.e. for r > R, satisfies 

g = S + gx + 0{r-P), V°g = 0{r~P-'), V°V'^5 = ^(r-^'-^) (2.16) 

where S = Sijdx^dx^ is the Euclidean metric on M™^'^, V*^ is the metric connection of 
go = S + gx, and p> 0. The AE case, the special case of AKK with A; = 0, is of relevance 
as it gives the correct asymptotics for static Lorentzian spacetimes with no non-extremal 
horizon, continued to Euclidean signature so that Euclidean time is non-compact. The 
AKK case with k = 1 and Xk = is of physical relevance since it gives the correct 
asymptotics for a static non-extremal black hole continued to Euclidean time where Xk 
is the Euclidean time circle. The AKK case with k > 1 may describe the Euclidean 
continuation of static Lorentzian spacetimes with Kaluza-Klein asymptotics, possibly with 
a non-extremal horizon in which case X^ contains an factor. We will now treat all these 
cases A; > simultaneously. 

The inverse metric satisfies g~^ = gQ^ + 0{r^^). Now introduce coordinates x^ for 
g hy x^ = (x^,y"') where are coordinates on Xk- It easily follows that the Christoffel 
symbols T'^p = 0(r^P^^) if at least one of the coordinate indices belong to the Cartesian 
ones x*, and that F^^ = T{gx)hc + 0(r~^~^). Thus assuming our reference metric g is also 
AKK with the same gx, we easily see 

^^ = 0(r-^'-i), ct, = 0{r-^P-^), (2.17) 

and thus (/>—)• at asymptotic infinity. We note that in fact we did not need the precise 
value of p and the argument works for any p > 0. 

If there are no other boundaries/asymptotic regions, a simple argument now allows 
us to use the maximum principle to rule out non-trivial solitons. Take Ai minus the 
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end E = (M™"'^ — Br) x Xk where r > R. Denote this manifold Yr, and note that it 
is compact with (m — l)-sphere boundary dYr = S"^~^. The AKK asymptotics imply 
(pldYr ^ Cr~^P~^ for some positive constant C. The maximum principle applied to Yr thus 
implies cp < Cr"^^"^ everywhere within Yr- Now, suppose there is a non-zero value of (p, 
say (j)Q > 0, at some point q in Ai. By taking a large enough r we can always arrange q to 
be in the interior of Yr and for Cr~'^'P~^ < cpQ. We therefore have a contradiction and thus 
we must have i;^ = everywhere. 

If there are other boundaries/asymptotic regions, this argument implies that if a non- 
trivial soliton exists on Ai, then the maximum of cj) must reside on/in a different bound- 
ary/asymptotic region. 

For completeness we record how to write AE metrics in the coordinate system (2.5). 
First write the Euclidean metric in polar coordinates 6 = dr^ + r^dO^. Then define p = 1/r 
so that yO — >• corresponds to r — )• cxd. The metric (and similarly the refererence metric) 
can then be written in the form 




(2.18) 



together with requiring the 'Dirichlet' conditions, a = 1, = and hij = Qij at p = 0. 
In fact, more precisely from the above we deduce a = 1 + 0{ff), oji = 0{pP~^'^) and 
hij = ilij + 0{ff'^'^). This is simply generalized to the AKK case. 

2.2.3 Asymptotically locally hyperbolic manifolds 

Next we consider the case where A = —^^^^^ < and {Ai,g) is a non-compact asymptoti- 
cally locally hyperbolic (ALH) manifold. We start with the standard conformal definition 
(see [24]) which requires the existence of: i) a Riemannian manifold with a boundary 
(M, g) such that A4 is diffeomorphic to M — dM; ii) a defining function z on M, such that 
z > and g = z^g in the interior, and z = and dz ^ Q on dM. For Einstein manifolds 
this is then of course supplemented by Ric(g) = Kg in M. Instead we will consider the 
Einstein-DeTurck soliton equation (2.1) such that the DeTurck vector ^ is constructed from 
a reference metric g in the same class. 

First define a vector field Z in M, by the properties that on dM it is normal (with 
respect to g) and g{Z, Z) = £^ on dM, such that Z^d^j^z = 1 in M. Let X be any vector 
field tangent to hypersurfaces of constant z (which includes dM), so one has X^d^z = 
and [X, Z] = 0. Since dz ^ on dM we may use z as a coordinate. We now define a notion 
of ALH, sufficiently general for the Einstein-DeTurck equation, by requiring the metric g 
satisfies 

g = i\dz'^ + h°) + 0{zP) (2.19) 

where is the boundary metric (an (m — l)-dimensional Riemannian metric), for some 
p > 0, together with derivative conditions 

(V|)'"i (V^)"^^ = 0(zP-"i) (2.20) 

where V" is the metric connection of g^ = i'^{dz'^ + h^) and are non-negative integers 
such that Til -\- n2 < 2. These conditions imply that the Riemann tensor oi g = z~'^g 
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approaches that of the standard hyperbohc space H'" as z — )• 0. More precisely, one can 
show that® 

I mc{g) -Ag\^ = Oiz^n \ Weyl(g) \' = 0{z'n (2.21) 

where the norms are of course taken with respect to the metric g. It is thus in this sense 
the metric is asymptotically hyperbolic. Observe that from the Ricci soliton equation it 
immediately follows that |V^p = 0{z^^). Shortly we will show that in fact the stronger 
statement = 0{z'^p) holds. 

Therefore, equivalently we may define our notion of an ALH manifold as having an end 
diffeomorphic to a quotient of — {z > zq}, such that g can be written for {) < z < zq as 
g = z~'^g with g given by (2.19) and (2.20). We will also choose the reference metric g to 
be in the same class with the same boundary metric. 

For practical purposes it is convenient to introduce an explicit general chart (z, x*) 
valid near the boundary, so Z = dz and X = di, which allows one to write g for small z > 
as 



a 



9 



{dz + ujidx")^ + hijdx^dx^ (2.22) 



where a,u}i,hij are all functions of {z,x^), z — )• is the asymptotic region, with — >• 
and hij — )■ h'^j. Given (2.20) it follows that the z — )■ behaviour of these functions is 



(2.23) 



{dzr{dkr{h,,-h%) = o{z^-'^^) 



where again rii are non- negative integers such that ni + 712 < 2. Note that in such coor- 
dinates we cannot assume smoothness (of g) near 2; = 0, since this property is not even 
satisfied by ALH Einstein metrics in the Feffer man- Graham expansion (they are differen- 
tiable to some finite order determined by m though). 

We remark that from the point of view of an elliptic boundary value problem we may 
regard taking this behaviour for a, uji and h as fixing Dirichlet data for all components of 
the metric g at the conformal boundary, i.e. a = = and hij = h^j at z = 0. 

As stated above we take the reference metric g also to be ALH in the above sense, with 
the same boundary metric in the conformal frame defined by z. Furthermore, we use 
the coordinate freedom in the chart {z,x^) and use Gaussian normal coordinates so a = ^ 
and Pi = 0. We thus just have 

hij = h% + 0{zP) . (2.24) 

The metric and reference metric may be recast to the form (2.5) by changing to the 
proper coord w = — logz and setting pi = Ui/z, jij = hij/z'^ and ^ij = hij/z"^. With our 
choice of asymptotic form for g and 5, it is then straightforward to check from (2.8), (2.9) 
that as z — 7- 

= 0{zn , l^'iri, = z'^h'^i^ij = 0{z'n , (2.25) 



*If we further assume that p > 2, as for C^-functions near the boundary, then one has | Ric{g) — Ag p = 
0{z^) and I Weyl(p) | = 0(2*), independently of p. 
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and therefore from (2.10) that 



cP = 0{z^P) . 



(2.26) 



It fohows that (/>—)• as z —7- 0. 

Now, assuming there are no other asymptotic regions of the manifold, and that there 
are no other boundaries, we can use an analogous argument to the one for the AE/AKK 
case above to rule out the existence of non-trivial Einstein-DeTurck solitons. Namely, 
consider the compact manifold with boundary given hy = M — E,^ where is the end 
corresponding to < z < e. The ALH conditions imply (j)\dXt < Ce^ for some positive 
constant C. The maximum principle applied to thus implies (j) < Ce^ everywhere in 
Xg. Now, suppose there is a non-zero value of (p at some point q in Ai. By taking small 
enough e we can always arrange q to be in the interior of and for Ce^ < (f>{q). This gives 
a contradiction and we thus must have (j) = everywhere. If there are other boundaries, 
these asymptotics imply a non-trivial Ricci soliton must have the maximum of (f> on another 
boundary or asymptotic region. 

2.2.4 Asymptotically extremal manifolds 

We are actually interested in static Lorentzian metrics. In this context, we are in particular 
interested in solutions with Killing horizons, such as black holes. In the Euclidean section, 
a non-extremal horizon simply appears as a smooth degeneration of the static Killing 
vector field d/dr (see §2.2.5), which is neither a boundary or an asymptotic end of the 
manifold. However, extremal horizons are qualitatively different, and from an Euclidean 
point of view should be thought of as asymptotic regions of M. More precisely, we will 
introduce the notion of asymptotically extremal manifolds, in the context non-compact, 
static Riemannian manifolds. These will be defined in such a way that their Lorentzian 
sections describe the exterior to a general static spacetime containing an extremal Killing 
horizon, as we now discuss. 

Consider a spacetime containing a smooth static extremal Killing horizon, of a Killing 
vector field V = d/dv. We will assume that the cross-sections of the horizon T-L are simply 
connected and that V is timelike just outside the horizon (as is the case for static black 
holes). It has been shown that the near-horizon geometry associated to such spacetimes 
can be written as [11] 



where ^) = is the horizon, ^'^j^ is the induced metric on T-L and To[x) is a positive function 
on Ti. The 2d metric in the square brackets is of course AdS2, and so the near- horizon 
geometry is simply a warped product of AdS2 and a Riemannian metric 7^^ on Ti . Although 
the above metric is strictly only valid for g > 0, g = is merely a coordinate singularity 
which may be removed by changing coordinates to {v, g) defined hy t = v + g~^ .'^ 

We are of course interested in the full spacetime metric near the extremal horizon, not 
just its near-horizon limit. It can be shown that outside the horizon we may write the 

^Standard Gaussian null coordinates are given by (i),r, a;") where r = To{x)q, see Appendix A. 




(2.27) 
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spacetime metric g as 

g = -T{q, x)g^dt^ + R{g, x) + guja{g, x)dx"-^ + -fabiQ, x)dx'^dx^ (2.28) 

for ^ > 0, where T,R > and T, R,uja,^ab are smooth functions of {g,x) at g = (at 
least), such that 

T|,=o = R\,=o = To{x) , jab\e=o = labi^) (2.29) 

, , , Ti (x) — i?i (x) 
X = ^ , ^ = constant 2.30 
Jo (a;) 

where Ti = 9gT|g=o and i?i = 9pi?|g=o. In Appendix A we prove this starting from a 
general smooth extremal Killing horizon written in Gaussian null coordinates, and then 
imposing staticity. 

It is worth noting that that given the metric (2.28), which appears singular at ^ = 0, 
it is easy to see the conditions (2.29) and (2.30) are sufficient for the existence of a smooth 
extremal Killing horizon at g = whose near-horizon limit coincides with (2.27). Explicitly, 
consider coordinate transformations of the form 

dt = dv+ ('^ + —]dg (2.31) 

for some constants ao,ai which are chosen to ensure that the metric functions and the 
inverse metric functions are smooth at g = 0. One finds that 

«o = -1 , ai = ^^^^ ■ (2.32) 

However, since oq, ai must be constants we learn a non-trivial condition that must be met 
to ensure smoothness of the extremal horizon, namely (2.30). Assuming this is the case 
the metric near g = then looks like 

g = To{x)[-g'^{l + 0{g))dv'^ + (1 + 0{g))dvdg] + IRuJadx^dg + -labdx^dx^ (2.33) 

which is indeed smooth and invertible at = 0. This allows one to extend to the region 
behind the horizon < 0. The surface ^) = is then a smooth extremal Killing horizon of 
d/dv as claimed. Finally, the near- horizon geometry defined by setting g ^ eg and t ^ t/e 
and taking e — )■ coincides with (2.27) as required. 

Next, consider Euclideanising the geometry by setting t = it. Now g ^ corresponds 
to a new asymptotic region. This motivates the following definition: we will say that a 
non-compact, static, Riemannian manifold {M, g^d / dr) is asymptotically extremal, if 
has an end in which the metric g can be written in coordinates (r, g, x°) such that 

g = T{g, x)g'^dT'^ + R{g, x) + gcoaig, x)dx''^ + "fab{Q, x)dx''dx^ (2.34) 

for sufficiently small > 0, where T,R > and T, R, oja, Jab are smooth at ^ = and the 
conditions (2.29) and (2.30) are satisfied. We may consider either non-compact or periodic 
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Euclidean time. We emphasize that in either case the Riemannian geometry will be smooth 
in the interior, even though in the periodic case the time circle degenerates as ^) — t- 0, and 
that by construction the Lorentzian continuation will be smooth at g = (and indeed one 
can continue through the horizon). 

Again, the class of metrics (2.34) may be cast in the form (2.5) by defining the proper 
coordinate w = log g and x* = (r, x"") and setting 

= R, f3idx' = gujadx" , 'jijdx'dx^ = Tg'^dr'^ + -iabdx°-dx^ . (2.35) 

We also assume the reference metric g is in the same class as g, with the same near-horizon 
geometry as g. That is, as g ^ Q the metric g can also be written as (2.34) for some T, R, 
^ab and (Da = (choosing normal coordinates), such that Tq = Rq = Tq and ^ab = lab- For 
such g and g, it is again straightforward to show from (2.8), (2.9) and (2.10) that as — )■ 

iu, = 0{g), er = 0, ia = 0{g), 4' = 0{g^) (2.36) 

and hence (/> — )• as ^) — )• 0. One can then use an analogous argument to the ALH case, 
to show the maximum principle rules out the existence of non-trivial solitons, unless there 
are other asymptotic regions or boundaries where (j) may achieve a maximum. 

2.2.5 Fictitious boundaries 

It is convenient to adapt coordinates to any isometries a geometry may have. Since there 
will be no dependence of the metric components on the coordinates corresponding to the 
symmetry directions, the effective dimension of the Einstein-DeTurck PDE system is re- 
duced. However, if the isometry does not have a free action, these adapted coordinates 
will not cover the fixed points. The canonical example of this is polar coordinates which 
strictly speaking do not cover the origin point. One possibility is to use other charts to 
cover the neighbourhood of these fixed points that are not adapted to the symmetry. How- 
ever, a more practical option is to treat these fixed point sets as 'fictitious boundaries' 
for the adapted coordinate chart, and to determine the necessary 'boundary conditions' in 
these adapted coordinates by deducing them from smoothness in a chart that does cover 
the fixed point set. Since these are not real boundaries of the manifold then of course our 
maximum principle implies that for a non-trivial soliton solution, the maximum of (p cannot 
reside at these fictitious boundaries. For a black hole the horizon and any axes of rota- 
tional symmetry are typical examples of submanifolds that are usefully treated as fictitious 
boundaries for coordinates adapted to static and rotational symmetry. We shall discuss 
these examples below, giving the fictitious boundary conditions to be imposed in these 
adapted coordinates, and also demonstrating how these conditions imply that a maximum 
of (j) cannot occur. 

Non-extremal horizons: First consider a static spacetime containing a generic smooth 
static non-extremal horizon, such that the static vector is timelike just outside the hori- 
zon. Outside the horizon, the metric may be written in coordinates adapted to the static 
isometry {t,'w,x°') and takes the form 

g = -T^w^dt"^ + W'^ {dw + w Qadx^f + labdx^'dx^ (2.37) 
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for u! > 0, where the horizon is at to = 0, such that T,W > and T, W, f^a, "jab are smooth 
functions of w'^ and x"" at w = and 



T 
W 



= K (2.38) 

w=Q 



where k > is a constant equal to the surface gravity. A proof of this is provided in 
Appendix B. Conversely, given the metric (2.37), which appears singular at w = 0, it 
is easy to see that the smoothness conditions just mentioned together with (2.38), are 
sufficient for the existence of a smooth non-extremal Killing horizon at w = with surface 
gravity k. To see this explicitly, consider the coordinate transformation 

r = w'^ , v = t+-logw (2.39) 

K 

in which case the metric near r = is 

g = — [-TKdv^ + dvdr] + 0{l)dr^ + W'^nadrdx" + {-fab + W^rna^b)dx''dx'' (2.40) 

with all functions smooth at r = 0. This shows that the metric at r = is smooth and 
invertible and allows one to extend to the region r < 0. Therefore as claimed the surface 
r = is a smooth non-extremal Killing horizon of d/dv with surface gravity k. 
Now, consider Euclideanising this geometry by setting t = it, so 

g = T^w^dT^ + {dw + w ^adx^f + -fabctx^dx^ (2.41) 

for w > Q such that T,W > and T, W, U,a, ^ab are smooth functions of w"^ and x"" aim = Q 
and (2.38) holds. The static Killing field d/dr now vanishes at w = Q where the isometry 
it generates has fixed action. However, due to the condition (2.38) this Riemannian metric 
is guaranteed to be smooth at w = 0, provided one identifies r periodically 

r ~ r + — . (2.42) 

This is easy to see by changing from the 2d "polar" coordinates {w,t), to Cartesian ones 
{x,y) defined by x = w cos kt and y = w sin kt, in the usual way. Therefore from the 
Euclidean point of view, a spacetime containing a static smooth non-extremal horizon 
corresponds to a smooth Riemannian manifold with a static U{1) isometry with a fixed 
action on a co-dimension-2 submanifold. Conversely, from the above it is easy to show 
that any smooth Riemannian manifold with a static U{1) isometry with a fixed action on 
a co-dimension-2 submanifold, has a Lorenzian continuation which can be extended to a 
static spacetime with a smooth non-extremal horizon. 

This metric may be cast in the form (2.5) by taking = (r, x") and setting 

a = W, 13a = w^a , -fijdx'dx^ = T^w'^dr^ + -fabdx^dx^ . (2.43) 

We take the reference metric g to be in the same class, so, 

a = W, %dx'dx^ = f'^w'^dT^ + ^abdx^dx^ . (2.44) 
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with the same condition that T, W, jab are smooth in x and atw = 0. From (2.8), (2.9) 
and (2.10) we see the smoothness of the functions T, W, ilaj 7a6 in w"^, and hkewise for the 
reference metric, imphes that = at it; = 0. This shows how the maximum principle 
imphes there can be no maximum of (p located at the horizon from the perspective of the 
"polar" coordinates. 



Axis of rotation: We note that although the U{1) above corresponds to Euclidean time, 
it could equally originate from a spatial rotational isometry. One can obtain an analogous 
result for a rotational SO{n) action with n > 2 that has fixed action on some submanifold. 
Then we may write the metric in "polar" coordinates as 

g = W'^ {dw + w Qadx^f + S^w^dn^ + jabdx^'dx'' (2.45) 

for w > 0, where dQ'^ is the unit (n — l)-sphere line element associated to the SO{n) action. 
Again, smoothness of the geometry implies W, S, Qq, jab are smooth functions of x and w'^ 
at w = 0, the only difference being that 

= 1 (2.46) 

since for n > 2 the (n — l)-sphere has curvature and must shrink at a particular rate. 
Again, taking the same behaviour for the reference metric, one finds d^cj) = at w = 0, 
and hence as we expect a maximum in (p is also excluded there. 



S_ 
W 



3. Ricci-DeTurck flow on manifolds with boundaries 

In this section we discuss Ricci-DeTurck flow as a method for solving the Einstein-DeTurck 
equation. 

Consider a one-parameter family of Riemannian metrics g labelled hy \, on M. The 
Ricci-DeTurck flow (with cosmological term) is a diffusive flow of the soliton equation (2.1), 
namely 

^ = -2 {R^, - V(^^,) - Ag^,) (3.1) 

where A is thought of as an auxiliary flow time. With the DeTurck choice of vector 
this flow is parabolic, and may be straightforwardly simulated numerically given "initial 
data". A natural choice we will use for this initial data is to take g\x=Q = g, where g is the 
reference metric used to define the DeTurck vector ^ (2.4), and hence ^|a=o = 0. If there 
are boundaries/asymptotic regions, then on top of the initial data, one needs to specify 
boundary/asymptotic conditions to make the problem well-posed, which we will discuss 
shortly. We note that it is well known that this flow is diffeomorphic to the Ricci fiow 
dg^y/d\ = —2 (^R^iu — Ag^Jj , where the diffeomorphism is generated by the vector field 
In the case of manifolds with boundaries, the metrics g and g will satify the same boundary 
conditions if ,^ = on the boundary (since ^ generates the diffeomorphism) . We note that 
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Ricci flows in the presence of certain Dirichlet boundaries have been previously considered 
in [4,5], and in particular in [5] a proof of uniqueness was given where symmetry reduced 
the problem to the flow of a cohomogeneity-one metric. For Neumann boundary conditions 
existence of Ricci flows has been proven [25] when the boundary is totally geodesic. 

Recall, that ultimately our motivation is to find new Einstein metrics. As pointed out 
in the introduction, there are two potential problems with using Ricci-DeTurck flow as a 
method to find Einstein metrics. The first is simply that if the solution we seek is not a 
stable fixed point of the flow, then in practice this method is harder to implement. The 
second is that if we do converge to a fixed point, then we in fact have a solution to the 
Einstein-DeTurck soliton equation. Therefore, if we are to use this flow as an algorithm to 
generate Einstein metrics, we would wish to avoid or at least identify non-trivial solitons. 
As we showed in the last section for a variety of boundary and asymptotic conditions on 
g of interest, a maximum principle forbids the existence of such Einstein-DeTurck solitons 
provided we also impose these conditions on the reference metric g (and indeed these 
boundary conditions all possess the property that ^ = on the boundary). 

The key issue in using Ricci-DeTurck flow in the presence of boundaries, is whether 
the boundary conditions required for the soliton equation are in fact preserved by the flow. 
This is of course a desirable feature since one can start with a metric which satisfies the 
boundary conditions of interest, and then if one flows to a fixed point one is guaranteed to 
have a solution satisfying the required boundary conditions. We will now discuss this for 
the types of boundary conditions introduced in 2.2. 

For general manifolds with a boundary we have discussed "modifled Dirichlet" and 
"mixed Neumann-Dirichlet" type boundary conditions in §2.2.1. In the modified Dirichlet 
case we have ^ = on the boundary. Since ^ generates infinitessimal diffeomorphisms along 
the Ricci-DeTurck flow, the vanishing of ^ at the boundary implies the boundary points are 
fixed during the fiow. Thus we expect this type of boundary condition to be preserved along 
the flow. In the mixed Neumann-Dirichlet case we again have that the normal component 
of vanishes, ensuring the boundary points are preserved by the flow. We note that in 
this case diffeomorphisms generated by are allowed along the flow that move the points 
within the boundary, since the tangential components of ^ are not restricted to vanish. 

Let us now address the case of flctitious boundaries considered in §2.2.5. In fact, 
the case of fictitious boundaries which arise from static non-extremal horizons, has been 
considered previously [1]. Indeed, it is easy to see the class of metrics (2.41) satisfying 
(2.38) are preserved by Ricci-DeTurck flow, using the fact these manifolds are equivalent 
to the class of smooth Riemannian manifolds with a static U{1) isometry with a fixed 
action on a codimension 2 submanifold. We simply note that since Ricci-DeTurck flow 
preserves isometrics, provided no singularities are met the U{1) action must leave the same 
codimension 2 submanifold fixed. An analogous argument works for the fictitious boundary 
that arises from an axis of symmetry. 

Now we consider the case where the manifold has an end in which the metric takes some 
asymptotic form. In fact the asymptotically Euclidean case has been considered before. In 
particular, Ricci flow on asymptotically Euclidean manifolds has been shown to preserve 
this class of manifolds [9]. Presumably this can be generalised to the Kaluza-Klein case. 
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We merely note here that since our sohton boundary conditions ensure Rf^i, = 0(r~^~^) 
and V(^^^) = 0(r~^~^) (see §2.2.2), as long as the (Ricci flat) metric gx on the internal 
manifold is independent of A, the Ricci-DeTurck flow equation is satisfied asymptotically.^*^ 
This leaves asymptotically locally hyperbolic manifolds and static asymptotically extremal 
manifolds. In the ALH case we will only provide a brief discussion, and note that re- 
cently, under assumptions of smoothness (which generally do not hold for even boundary 
dimension), Ricci flow has been shown to preserve ALH metrics [10]. The asymptotically 
extremal case has not been previously considered and the majority of the remainder of this 
section will be devoted to proving the flow preserves the class of asymptotically extremal 
manifolds. 

3.1 Asymptotically locally hyperbolic manifolds 

An important question which we now address is whether Ricci/Ricci-DeTurck flow preserves 
the class of ALH manifolds in question. This is crucial in order to guarantee that given an 
initial ALH manifold, any flxed point of the flow is also an ALH manifold. Bahuaud has 
recently proved short time existence of Ricci flow that preserves smoothness [10].^^ The 
metric can be written as 

g = ^ {dz^ + {hij{z,x) + zvij{X, z,x)) dx^dx^^ (3-2) 

for a flow time A, so hij is fixed and Vij varies such that Vij\x=Q = 0. Hence hij gives the 
conformal boundary metric which is flxed in flow time. Then Bahuaud has shown that 
provided z'^g is smooth up to the boundary initially, then the Ricci flow exists for short 
times, and z'^g remains smooth up to the boundary. 

However, in even boundary dimensions the Feffer man- Graham expansion shows that 
an Einstein metric need not be smooth. Nevertheless, it seems reasonable to expect that 
even in the non-smooth case, provided the boundary metric is fixed, Ricci/Ricci-DeTurck 
flow does indeed preserve the class of ALH manifolds as defined above in section 2.2.3. 
To be more precise we will define ALH Ricci-DeTurck flow in an analogous manner to 
the Einstein-DeTurck equation in §2.2.3. That is supplement properties (i) and (ii) of the 
standard conformal definition of ALH at the beginning of §2.2.3, together with the defining 
function z being independent of flow time A, and of course imposing the Ricci-DeTurck 
flow (3.1) with A < on A^. The key question is then whether the remaining boundary 
conditions given in §2.2.3 are preserved by the flow. Assume the boundary conditions in 
§2.2.3 are satisfied along the flow, so asymptotically as z — )■ we have I{ic{g) ~ Ag and 

— )• 0, and therefore from the flow equation dg/dX — )■ 0. Therefore we learn that a 
necessary condition that the flow is well defined is that the conformal boundary metric 
is fixed along the flow, i.e. it is independent of A. It would be very interesting to show this 
is also sufficient in the case that the metric is not smooth. 

An interesting question is whether there is some analog of the Fefferman-Graham 
expansion that is valid along the fiow, and possibly some notion of boundary stress tensor. 

^"in fact naively this argument suggests that the ADM mass does not change along the flow, so one has 
to take these arguments with caution, see [9]. 

^^We thank Eric Woolgar for bringing this reference to our attention. 
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In the Appendix C we argue that a naive application of the Feffer man- Graham expansion 
can not apply along the flow, by studying the simpler problem of diffusion of a scalar in 
hyperbolic space. It would be very interesting to investigate this further, and determine 
whether there are suitable asymptotic expansions and boundary stress tensors. 

3.2 Static asymptotically extremal manifolds 

Now we address whether Ricci-DeTurck flow preserves the static asymptotically extremal 
manifolds, which we defined in §2.2.4. This question has two parts. The first is simply 
whether the flow in the class of metrics (2.34) preserves the conditions (2.29); if so, the 
second question is whether the smoothness condition (2.30) is preserved by the flow. 

Let us address the first part of the question, i.e whether the flow preserves the near- 
horizon geometry. Thus, suppose (2.34) now depends on A, such that at A = the near- 
horizon geometry coincides with (2.27). Explictly, the "near-horizon" limit is extracted by 
setting Q ^ eg and r — >• r/e and e — )• 0, so that the full metric (2.34) becomes: 

go = Q^Toix- \)dT^ + ^"^^'^^"^^^ + 7°b(x; \)dx-dx^ (3.3) 

subject to the "initial" conditions To(x;0) = Tq{x) = Tq{x), Rq{x]Q) = Tq{x) = Tq{x) and 
lab{x; 0) = 7°fe(x) = 7°b(a;), i-e- 9o|a=o = Qnh = go- The near-horizon limit of the fuU flow 
equation gives 

dX 

where .^o is the DeTurck vector associated to the metric go and background metric 
(which is the near-horizon limit of the full DeTurck vector). This can be thought of as 
a Ricci-DeTurck flow equation for the near-horizon geometry itself. Now, let us assume 
that our initial near- horizon geometry solves the Einstein equations (Ric((7o) — ^go)\\=o = 
0. Also note that since the initial near- horizon geometry coincides with the reference 
one we must have ^o|a=o = 0. Therefore, it is clear that a solution to the near-horizon 
geometry flow equation (3.4) for A > is given by taking To(x; A) = To{x) = Ro{x; A) and 
^abi^'i ~ 7afe(^)- Now, appealing to uniqueness of Ricci-DeTurck flow,^^ this must be the 
only solution with such initial conditions. This shows that if one starts with a near-horizon 
geometry which solves Einstein's equations, then it will not change under the flow. In other 
words the flow preserves the near-horizon geometry. From now on we will assume this to 
be the case. 

The second question which needs to be addressed is whether the smoothness condition 
(2.30) is preserved along the flow. This is a more non-trivial question which requires 
expanding the full Ricci-DeTurck flow equation to first order in g, to derive flow equations 
for Ti, and hence ip deflned above. The details of this calculation are given in Appendix 
A. We flnd a remarkably simple flow equation for iIj onT-L: 

dih „9 dTn , , , , 



-2 (mcto) - \Ci„go - Aso) (3.4) 



^^This has been shown for certain complete non-compact Riemannian manifolds [26,27] and we assume 
it to be the case here. 
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where V and • are the Levi-Civita connection and contraction with respect to the metric 
7^^. Recah that smoothness of the extremal horizon requires ■0 to be a constant (2.30). 
If we take an initial condition V'|a=o = V'o where i/jq is a constant, then well-posedness 
of the flow equation for tp shows that ip = iI:q for all A > 0. This is simply because, 
ip = il^Q clearly solves the flow equation with such an initial condition, and by uniqueness 
this must be the only such solution. This proves that Ricci-DeTurck flow preserves the 
class of static asymptotically extremal Riemannian manifolds, or in Lorenzian language, 
the class of smooth extremal static horizons with a given near-horizon solution to Einstein's 
equation. 

4. Example: AdSs gravity dual to strongly coupled CFT4 on a Schwarzschild 
background in the Unruh or Boulware vacua 

In this section we will demonstrate the use of the Ricci-DeTurck flow to find a f)d Einstein 
metric with A < 0. This will be a static axisymmetric spacetime, with asymptotic AdS 
boundary metric being Ad Schwarzschild, and with metric asymptoting to the AdS Poincare 
horizon far from the boundary. In the context of AdS/CFT, taking a trivial product with 
a round then gives the classical bulk metric dual to A/" = 4 super Yang-Mills on a 
Schwarzschild background. The IR boundary conditions we impose (i.e. that the geometry 
tends to the AdS Poincare horizon) are consistent with the classical geometry describing 
the dominant bulk saddle point dual to the the GET in the Unruh or Boulware vacua, 
where asymptotically the theory is in vacuum. We note that they are different to those 
considered in [13], who conjecture the behaviour of the dual geometry for the Hartle- 
Hawking vacuum, where one has a finite temperature horizon in the IR at equilibrium with 
the boundary horizon, and one expects the GET to asymptotically be at finite temperature. 
Here we have chosen to study the classical dual describing the GET Unruh and Boulware 
vacua as, i) it nicely illustrates the implementation of boundary conditions for extremal 
horizons discussed formally above, and ii) the physics of the Unruh and Boulware vacua is 
remarkably different for this strongly coupled GET than for free field theory. This latter 
point has important implications for the related issue of existence of Randall-Sundrum 
braneworld black holes [13,28-33]. 

As discussed above we treat the static solution by considering the Euclidean section, 
which is a Riemannian manifold with two asymptotic regions (one locally hyperbolic, the 
other extremal) and smooth in the interior. We will find that a single coordinate chart 
suffices to cover this manifold, and in the following sections will discuss and motivate its 
choice. As we shall see, the problem is static and axisymmetric, and hence of cohomogeniety 
two. It depends non-trivially on two coordinates, essentially the boundary (GET) radial 
coordinate, and the bulk coordinate, and compactifying these the problem is formulated on 
a rectangular domain. Two boundaries are fictitious, associated to the axis of rotational 
symmetry and the vanishing of d/dr at the horizon. The remaining two are the asymptotic 
AdS boundary (corresponding to the UV of the GET) whose conformal boundary metric 
we specify as Schwarzschild, and the asymptotic extremal horizon (corresponding to the IR 
of the GET) we choose to have a near horizon geometry given by Poincare AdSs. Using our 
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discussion above we are able to rule out a maximum of (f) on all of these boundaries, and 
hence the existence of any non-trivial Ricci solitons. Crucially, this implies that provided 
the DeTurck flow converges to a fixed point, it must give an Einstein metric relevant for 
our problem. Finally, we numerically simulate the flow and show that indeed a stable flxed 
point exists. 

4.1 Static and axisymmetric metrics: coordinates for the problem 

First let us consider AdS5 in Poincarc coordinates 



9AdS, = ^ [dz^ + dR^ + R^d^lr^ - dt") (4.1) 



where z > and i? > (equality is not strictly allowed in these coordinates, and corre- 
sponds to the the conformal boundary and the fixed point set of 50(3) respectively), and 
A = — ^. Now introduce new coordinates (r, re) by 



with < X < 1 and < r < 1 to ensure an invertible transformation. The metric in these 
coordinates is 

£2 / 47-2 4 \ 

where for covenience we define 

/(r) = 1 - , g{x) = 2 - x\ (4.3) 

In these coordinates the conformal boundary is in the asymptotic region x — >■ 1, and .x = 
is the fixed point set of the S0{2>) symmetry. Moreover, these coordinates are particularly 
adapted to the Poincare horizon which is at r = 1 and x < 1. Indeed, for x < 1, we can 
define new coordinates (w, q) by /(r) = q and t = v + ^ to get 



9AdS5 



{1-x 



2\2 



-Q dv +2dvdg-\ — —- dx +xg[x)d^(n\ 
g{x) ^ \ 



(4.4) 



which confirms that = is a smooth extremal Killing horizon w.r.t. d/dv, and allows 
one to extend through the horizon to the region ^ < 0. This expresses AdS^ as a warped 
product of AdS2 and a metric jab on cross-sections of the Poincare horizon, i.e. as a static 
near-horizon geometry with non-compact cross-sections of the horizon. The geometry 
induced on the Poincare horizon is thus given by 

^2/4 \ 
7„6dx"dx^ = dx^ + dnf^^j (4.5) 



^^Note that the general form for a static near-horizon geometry was derived in [11], and we have just 
written Poincare AdS in this general form. 
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where < x < 1, which is easily seen to be locally conformal to the standard round metric 
on (by changing coordinate to y = 1 — x^). Hence, topologically one can think of 
spatial sections of the Poincare horizon as open hemispheres of S^. Note that as x — >• 
the metric is of course regular and looks like the origin of M'^ , whereas x — >• 1 corresponds 
to an asymptotic region (which touches the conformal boundary). 

We wish to find a static solution which is asymptotically locally AdS, with the metric 
on the conformal boundary in the same conformal class as 4d Schwarzschild, such that it 
approaches the AdS Poincare horizon far from the conformal boundary. In the Euclidean 
section our manifold is thus ALH with conformal boundary metric given by Euclidean 
Schwarzschild which of course has U{1) x SO{3) isometry. In fact, it has been shown that 
regular static ALH Einstein metrics must inherit the isometry of the boundary metric [34] . 
Therefore our full spacetime metric must have U{1) x 50(3) symmetry. 

We now write an ansatz for such a static Riemannian metric in the coordinates devel- 
oped above 

a. - ^ (r'Tdr' + ^^^dn' + dr^ + dx^ + ^ " drdx\ 



z{r,x) = ^ — ^ (4.6) 
1 — 



where f,g are defined by (4.3), X = {T, S, A, B, F} are smooth functions only depending 
on the coordinates (r, x), thus making the static isometry generated by d/dr and the 
spherical symmetry both manifest. As above, the coordinate ranges are < x < 1 and 
< r < 1. Eventually we will treat the asymptotic regions x = 1 (conformal infinity) and 
r = 1 (Poincare horizon) as boundaries and work on the square 0<x<l,0<r<l 
with X being bounded on this domain (and smooth everywhere except at x = 1). The 
other domain boundaries r = and x = we will consider as fictitious ones associated to 
the (non-extremal) black hole horizon and the axis of rotational symmetry. Of course, 
we require that T,S > and AB > r'^x^g{x)F'^ / 16 (i.e. the determinant of the r, x part 
of the metric is positive), to ensure we have a good Riemannian metric. We note that for 
r^T = S = A = B = 1 and F = the metric reduces to AdSs in Poincare coordinates as 
written in equation (4.3) ~ however in this case T is not smooth at r = and hence AdSs 
is not strictly within our ansatz. 



4.2 Boundary and asymptotic conditions 

Having defined our coordinate domain, we now discuss the behaviour we impose on the 
metric functions at its various boundaries. Each boundary has been discussed in a general 
setting in the earlier §2.2, and we now apply those general results to our specific problem. 

^*We note that after choosing the reference metric there is in principle sufficient coordinate freedom in 
the Einstein solutions of the Einstein-DeTurck equation to choose these boundaries of the coordinates [1]. 
This is reflected in the Ricci-DeTurck flow by the vanishing of the normal component of ^ at the edge of 
this domain which follows from the various boundary conditions. 
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A defining feature of our metric is that it should be ALH with the boundary metric 
conformal to Euclidean Schwarzschild 

hscku. = {^-^) + + R'dnf,^ (4.7) 

where r ~ T+47ri?o is required for smoothness. This can be written in the above coordinates 
by setting R = Ro/{l — r'^) (so r = is the horizon) which gives 



hschw 



(4.8) 



f{ry 

where /(r) = 1 — as above. It is clear from our ansatz (4.6) that if we impose 

rI, B^f, T^l, RI, F^O, as x^l (4.9) 

then 

which indeed shows that the conformal boundary metric is conformal to the Schwarzschild 
metric as required. Recall we do not assume smoothness of X at x = 1, rather the 
asymptotic conditions discussed in (2.2.3). 

Following our general discussion of extremal horizons in §2.2.4, since we wish our metric 
to approach the Poincare horizon (4.3) smoothly, we must impose 

A = e^{\- r)Ai , B = e^{\- r)Bi , (4.11) 

T = f + {l-r)Ti, S = i^ + {1- r)Si , F = (1 - r)Fi , 

Ti\g=Q — Ai\g=Q = constant . (4-12) 

where Xi = {Ti, Si, Ai, Bi, Fi} are all smooth at r = 1 and a; < 1, i.e. the Poincare 
horizon. As shown in §2.2.4, these conditions are necessary and sufficient for our spacetime 
to contain a smooth extremal Killing future horizon at r = 1, a; < 1 which coincides 
with the standard Poincare horizon of AdSs. To see this explictly, one must pass to the 
Lorenzian section t = ir and consider the coordinate transformation defined by = 1 — 
and (2.31) and (2.32). Indeed, if the boundary conditions (4.11) and (4.12) are met the 
metric near ^ = takes the form 

f 

= [-Q^{1 + 0{Q))dv^ + 2(1 + 0{g))dvdQ + xO{l)dQdx 



(1 -X2)2 



4(1 + 0(g)) dx"^ , , 2/-, ^/ ,r^2 

+- 7^ — + g{x)x^{i + 0{Q))d^%. 

g[x) 



(4.13) 



thus showing = is a smooth extremal Killing horizon of d/dv which coincides with that 
of the Poincare horizon. 

It is worth mentioning at this stage that in fact (4.4) is not the most general boundary 
condition one could take at a static extremal horizon. As is well known, the near-horizon 
geometry of extremal Killing horizons is in fact completely determined by the Einstein 
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equations (see e.g. [35]). Therefore, to determine the most general boundary conditions at 
such a horizon in this situation one would have to provide a classification of 5d vacuum 
static near-horizon geometries with an SO{3) symmetry. This problem has been considered 
in [36] . It is worth noting that static 4d vacuum near- horizon geometries can be completely 
classified (with no assumptions of symmetry) [37] and turn out to be a 1-parameter gen- 
eralisation of AdS4 written in Poincare horizon adapted coordinates. The results of [36] 
indicate a 1-parameter family also in 5d. Thus, although more general static near-horizon 
geometries in our 5d case are likely to exist, since we are interested in the application to 
AdS/CFT, we restrict our attention only to the simplest case of the Poincare AdSs near 
horizon geometry. 

4.3 Ricci-DeTurck flow and the absence of solitons 

Our Riemannian manifold has two ends and it is ALH in one end (with conformal boundary 
metric in the Schwarzschild class) and asymptotically has a Lorentzian continuation to the 
extremal Poincare horizon of AdS in the other end. We must also ensure that the metric 
is everywhere smooth in the interior of the manifold. Two regions where this can fail are 
where the U{1) and 5*0(3) symmetries have fixed points. These correspond to r = the 
non-extremal horizon arising from the Schwarzschild metric on the conformal boundary, 
and X = the axis of the spherical symmetry, respectively. To discuss these regions 
properly we introduce fictitious (from the full manifold point of view) boundaries at r = 
and X = and work on the square {0<r<l, 0<x<l}. Topologically, this square can 
be identified with the orbit space O = M./[SO{3) x U{1)], which has the natural structure 
of a manifold with boundaries. 

As discussed in §2.2.5, at the horizon r = smoothness of the manifold requires the 
metric functions T, A, B, S, F to be smooth in near r = 0, and T = k^A at r = 0, where 
K is a constant so that the period of r is 27r/K. The Schwarzschild conformal boundary 
condition given above in (4.9) determines k = 1/{2Rq). Likewise, at the axis of spherical 
symmetry, we require T, A, B, S, F to be smooth in near x = 0, and in addition, S/B = 1 
at X = 0. 

We consider Ricci-DeTurck flow in this class of manifolds, with g given as in (4.6) 
in terms of the functions T,A,B,S,F with behaviours described above. We choose the 
reference metric g to be in the same class, and in normal coordinates to all boundaries, as 
we chose in the earlier §2.2. We make the specific choice 



9 = 9 



(4.14) 

T,A,B,S=1,F=0 



Whilst we have proved that the metric behaviour near the fictitious boundaries and the 
extremal horizon is preserved under the Ricci-DeTurck flow, we have no formal proof that 
this is so for the ALH, although it seems reasonable to expect this. We emphasize that 
it would be interesting to try to prove that this is so, and merely observe later in our 
simulations that it appears to be the case. 

Suppose we flow from some initial metric in our class to a fixed point in it. We have 
argued in §2.2 that no maximum of (j) can occur in the ALH and extremal asymptotic regions 
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where (p ^ 0, and since the other boundaries are fictitious no maximum can occur there 
either (as discussed, the normal derivative of (p vanishes on these fictitious boundaries). 
Then our earher arguments applied here demonstrate there can be no non-trivial solitons. 
This is a powerful result. Whilst one might have worried that many solitons could in 
principle exist and would then complicate using Ricci-DeTurck flow to flnd the metric of 
interest, we have the beautifully simple result that in fact there can be no solitons at all. 
Hence we conclude that in our case any fixed point of Ricci-DeTurck flow must be an 
Einstein metric, and therefore precisely relevant to the AdS/CFT problem of interest. 

4.4 Numerical simulation of Ricci-DeTurck flow 

There are two remaining issues that could complicate using DeTurck flow to generate the 
Einstein metric we are interested in. Firstly, we have no proof that an Einstein metric 
satisfying the various conditions exists. Since we have ruled out the existence of solitons, 
this would presumably have to manifest itself in the DeTurck flow inevitably reaching a 
finite time singularity for any initial data. Secondly, black hole fixed points have previously 
been shown to often be unstable under Ricci flow [38] , due to the famous Gross- Perry- Yaffe 
Euclidean negative mode and its generalizations [12]. However, in our case here the black 
hole metric is in part flxed at the AdS conformal boundary, and one then hopes that this 
projects out any putative negative modes. We have no way to address these questions 
except to simulate the DeTurck flow and observe what happens. As we shall see, we find 
an elegant and simple flow that quickly asymptotes to a fixed point, obviously indicating 
its stability and existence. 

There are two scales in the problem, namely A and Rq. We choose units so that the 
AdS length scale £, defined as A = —4/^^, is one. Since the remaining scale Rq enters 
only in the conformal boundary metric, and we may always rescale this since it is only the 
conformal class that is relevant, we may without loss of generality choose this to be one in 
our units. Having fixed £ = 1 and Rq = 1, where is no remaining data in the problem. 

Given that we have a reference metric g, a natural starting point for the flow is to take 
initial data for the flow g{X = 0) = g. All results presented later are for this choice. We 
have tried other choices, and flnd similar qualitative behaviour of the flow, and precisely 
the same fixed point. 

We truncate the continuum DeTurck equations by specifying the metric functions 
T, A, B, S, F on a finite set of points in the coordinate domain 0<r<l,0<x<l. We 
must then implement the various behaviours on the boundaries of our domain. Firstly at 
the AdS boundary x = 1 we simply impose the Dirichlet condition T = A = B = S = 1 and 
F = 0. Note our initial metric and reference metric, both g given in (4.14), satisfy the ALH 
boundary (and fall off) conditions we defined in §2.2.3. Thus from the discussion in §3.1, we 
expect this ALH behaviour will be preserved along the fiow. Indeed this is what we observe 
in the simulations. Secondly at the extremal Poincare horizon, r = 1 and x < 1, we impose 

^^This can be seen explicitly by defining a dimensionless angular coordinate 6 = r/(2i?o) (which has 
period 2tt) in terms of which hschw has an overall factor of Rq. Then in equation (4.10) one gets a 
dimensionless factor of Ro/£^ in the conformal factor of the boundary metric. By then rescaling 1 — a; it is 
always possible to set Ro/£ = 1. 
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the Dirichlet condition T = A = B = S = 1 and F = 0. As discussed earlier, our choice 
of g and hence initial and reference metric asymptote smoothly to these values in (1 — r) 
and Ti — is a constant. Furthermore, since we proved these properties are preserved 
by the flow for a general static extremal horizon whose near-horizon geometry solves the 
Einstein equation, we may apply this to the case at hand and deduce the flow preserves 
the near-horizon geometry and smoothness of the corresponding extremal horizon. 

At the fictitious horizon boundary, r = 0, the smoothness in implies we may simply 
determine the metric functions at the boundary points from their values at interior points 
by imposing Neumann boundary conditions on the metric functions, drT = . . . = drF = 
at r = 0. We note that the initial data and reference metric are indeed smooth in r^, 
and satisfy the condition T = k^A at r = 0, where k = l/(2iio) = 1/2 for our choice of 
parameters. Since the DeTurck flow will preserve this smoothness in and the condition 
T = K?A, in principle there is nothing more to do. However, we have found that directly 
imposing T = k?A improves numerical accuracy and stability, essentially as it eliminates 
potential 1/r divergences that arise from discretization error. 

The remaining domain boundary at x = is the fictitious axis of symmetry boundary. 
The required behaviour is smoothness of metric functions in x^, and S = B aX x = ^. This 
may simply be used to determine the values of the metric functions at boundary points from 
their values in the interior from Neumann boundary conditions, dxT = . . . = dxF = at 
X = Q. Again we observe that since is smooth in x'^ and obeys S = B ai x = the initial 
metric and reference metric have the desired regular behaviour for the axis of symmetry, 
and DeTurck flow will preserve this. As with the horizon we have found that directly 
imposing S" = on the axis significantly improves numerical accuracy and stability. 

We simulate the flow using two independent codes, both based on forward Euler dif- 
ferencing in flow time; one code uses a second order flnite difference in r, x (up to 160 x 160 
points), and the other a quasi-spectral representation of r, x (up to 40 x 40 points). We in- 
deed flnd a stable fixed point for the numerical flows, and the metric functions T, A, B, S, F 
remain smooth in the interior of our coordinate domain, and compatible with our boundary 
conditions. In Fig. 1 we plot the maximum value of the flow time derivative, d/dX, of 
one of the metric functions, T, over our coordinate domain. We see that the function T 
approaches the fixed point as a power law which we estimate as ~ with p ~ 2. We see 
analogous scaling towards the fixed point for the remaining metric functions, and all other 
geometric quantities we have computed. 

We are ultimately interested in the metric at the fixed point. Using Ricci-DeTurck 
flow to find this there are two sources of numerical error. Firstly, and most fundamentally, 
having spatially discretized the solution, the fixed point will solve an approximation to the 
Einstein-DeTurck equation, not the continuum equation. Secondly, if one simulates the 
flow for a finite time one will not actually reach the fixed point. The first source of error 
is inescapable, and is fundamental to any numerical approach. However, we can avoid 
the second source of error by making sure that we fiow for long enough so that the error 
introduced by terminating the fiow at finite flow time is small compared to the actual 
discretization error. 

With our boundary conditions, the maximum principle guarantees that (j) should vanish 
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at the fixed point of tfie flow in tlie continuum, and therefore we can use this quantity to 
estimate the numerical discretization error for the flow. In Fig. 2 we show the evolution of 
the maximum value of (j), denoted by 4>ma^, as a function of flow time. This is initially zero 
since we are showing results for a flow where the initial data was chosen to be the reference 
metric g. It evolves to be non-zero, although interestingly it is never very large. At late 
times as W6 a.pproa>ch the fixed point in tlie continuuni W6 would expect (/'max 

— (since 

the fixed point must have ^ = 0). However, the numerical discretization error implies that 
at the numerical fixed point, the scalar <j) will be non-zero. Taking higher resolution, the 
non-zero value of (pmax at the numerical fixed point should go to zero, its continuum value, 
scaling in accord with the nature of the discretization. The figure shows this behaviour 
quite clearly. Around A ~ 1.5 the value of (/)max saturates for the 20 x 20 quasi-spectral 
flow, indicating that whilst the metric functions are still evolving as we are at finite fiow 
time, the error in the solution is now dominated by the numerical discretization error. 
Increasing to 30 x 30 and then 40 x 40 resolutions we see that one must fiow for longer 
before the finite flow error becomes subdominant to the discretization error. In Fig 3 we 
plot the saturating value of (prnax as a function of resolution for the quasi-spectral flows, and 
find beautiful exponential convergence in the number of grid points as we should expect 
provided T, A, B, S, F are suitably smooth functions of the coordinates. 

The square root of 0max provides an indicator of the global fractional numerical error 
(since (j) is contructed from the square of ^) , and we see that for quasi-spectral 40 x 40 this 
error is already impressively small, < 10~^, and much smaller than our second order finite 
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Figure 2: Maximum value of 4> in the whole domain as a function of the Ricci flow time A for 
different spatial resolutions obtained using the quasi-spectral code. The blue curve corresponds to 
the 20 X 20 data, the red one to the 30 x 30 data and the black curve to the 40 x 40 data. The value 
of (/>max at the fixed point (constant section of the curves above) decreases as the spatial resolution 
is increased. 



difference code achieves at the highest resolution implemented. Data presented from this 
point on will be for the quasi-spectral solution obtained by flowing for sufficient time that 
the numerical discretization error dominates. 

In fact we have also used the Newton method, discussed in [1], to directly solve the 
Einstein-DeTurck equation to find the fixed point using identical spatial discretizations. We 
find solutions using both second order and quasi-spectral codes that are precisely consistent 
with those found by simulating the Ricci-DeTurck flow for long times. We emphasize that 
implementing the Newton method is considerably more complicated than implementing 
the flow method, although the former is considerably quicker taking hours rather than 
days to compute the 40 x 40 quasi-spectral solution running in Mathematica on a single 
processor machine. However we should highlight that the path taken through the space of 
geometries by the Ricci-DeTurck flow method is interesting, and as emphasized in [1] does 
not depend on the choice of reference metric. This should be contrasted with the Newton 
method where the iterations that improve the solution inherently depend on the choice of 
reference metric and therefore contain no geometric information. 

We may embed the horizon of our geometry along the flow into hyperbolic space with 
the same AdS length i, in axisymmetric coordinates, ie. ds'^ = || ^dz^ -|- dR^ + R'^di}'^^)^^ 
as the surface of revolution R = R{z), such that the induced metric of this surface is the 
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Figure 3: Saturating value of </>max as a function of the number of grid points A'' (in either the r or 
X directions) for the quasi-spectral code. 0max converges to zero exponentially with N, as expected 
for a suitably smooth solution. 




z 



Figure 4: Embedding into hyperbolic space, ds^ = ^ ydz^ + dE? + R^dfl^2)j ; of the spatial cross 
sections of the horizon along the flow as curves R{z). The dashed line corresponds to the initial 
data, for which the horizon is round, and the thick black line is the embedding of the horizon of 
the fixed point. The snapshots are drawn at intervals of A of 0.05. 

same as that of the horizon. The freedom in the embedding is fixed by setting R{0) = 1. 
Fig. 4 then shows the evolution of this surface from the initial data to late times in the 
flow. We note that the maximum extent in z of the embedding approaches its fixed point 
value in the same manner as the metric functions discussed above, namely as ~ X~^. 
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We see that the shape of the horizon is very different from that of a black string, 
ds^ = 1^ [dz^ + hschw), where hschw is the Ad Schwarzschild metric (4.8). This is also an 
Einstein metric, but does not have a regular Cauchy horizon, and is unstable to a Gregory- 
Laflamme instability [39,40]. From our embeddings we see the horizon of our solution does 
not look string-like for z > 1 {a string would have embedding R{z) = 1). Following the 
discussion in [31] , one would expect a Gregory-Laflamme instability only if the embedding 
behaved as R{z) ~ 1 up to quite large z, say z > 5 before 'capping off', otherwise a 
wavelength of the instability cannot fit into the string- like portion (see Fig. 2 in [31]). 
Since this is not the case, we do not expect our solution is unstable to a Gregory-Laflamme 
instability. We emphasize here that since our flow stably approaches the fixed point, our 
solution cannot have negative modes which respect the static axisymmetry. In fact we 
have explicitly confirmed this by a direct calculation of the spectrum of the linearization of 
the Einstein-DeTurck equation about the solution assuming its symmetries. Whilst there 
is no direct correlation between Euclidean negative modes and Lorentzian instability, it 
would have been suggestive of a potential instability had negative mode(s) existed. Thus 
we find nothing to suggest dynamical instability, but emphasize that one must perform 
a Lorentzian linear analysis to actually determine dynamical stability. We have not done 
this, and it would obviously be interesting to do so. 

Finally, in Fig. 5 we depict the evolution of the curvature invariant {C^upaC^'^''")^l 
along the fiow. As this figure shows, all along the fiow the Weyl tensor vanishes at both the 
Poincare horizon and at the AdS boundary. This is consistent with the Ricci-DeTurck fiow 
preserving ALH metrics with a fixed boundary metric, as we discussed earlier in §3.1. At 
the fixed point, from the behaviour of this curvature invariant near these two asymptotic 
regions of the spacetime we can estimate the behaviour of the Weyl tensor, and we find that 
Cpupa ~ (1 — near the boundary of AdS and C^^pa ~ i^ — f) near the Poincare horizon. 
We have also checked that for our solution the condition Ti\r=i — ^i|r=i = constant is 
satisfied (within numerical accuracy). Recall from the discussion around eq. (4.12) that 
this is a necessary condition for the regularity of a static extremal Killing horizon (given 
the near- horizon geometry is). For the initial data (and hence our reference metric g) the 
curvature is quite small everywhere, but along the fiow it grows near the horizon and the 
axis of symmetry. It is worth noting that for the fixed point solution the curvature is finite 
everywhere inside our coordinate domain. 



4.5 Boundary stress tensor 



Now let us consider the boundary stress tensor. For a fixed point with = 0, one may 
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Figure 5: Evolution of (C^i^po-C'"''"^)!^ along the flow. The snapshots correspond to A = 
0,0.1,0.2,0.5,1 and the fixed point. As shown in the plot, the Weyl tensor remains zero at the 
Poincare horizon (r = 1) and also at the AdS boundary {x ~ 1). The latter suggests that the ALH 
boundary conditions are preserved by the flow. 



compute the asymptotic behaviour of the functions X, and find, 



T = 1 - (1 - x'Yil - r') + (1 - x'f U{r) + ..., 
8 = 1 + ^(1- - + (1 - xY S4(r) + . . . , 

A = l- {1- x^f{l - r^) + ^(1 - x^)'^[3 - Sr^ + 5r^ - 4t4(r) -8s4{r)] + ... , 
i? = l + i(l-xY(l-r4) + ... , 

F=-{1- x^il - r2) + ^(1 - r2)2(log(l - x^){l - x^f) + Mr){l - x^ f + . . . , 

(4.15) 

where t4(r) is data to be determined by fitting from the bulk solution, and s^^r) is given 
in terms of this from, 



8(1 - 3r2)(s4 + ti) - (1 - r^) [3 - 27 + 30 - 4 r(2 ^ + t'^}] = . 



(4.16) 
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After a coordinate transform to the usual Feffer man- Graham form 

7i 
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dz'^ + hschw + z^hi + 0(z^ 



(4.17) 



with Schwarzschild conformal boundary metric hschw given by (4.8), one can extract the 
leading large Nc behaviour of the vev of the boundary stress tensor, {T--'), 



1 1 /3i2o , . 



{U{R) + 2 S4iR)) 



<R 



8R 



(4.18) 



where the first two terms correspond to the rr and RR components respectively. Equation 
(4.16) implies the stress tensor is conserved. Notice that the stress tensor is traceless, as it 
should be since the boundary metric is Ricci flat and therefore there is no anomaly in the 
boundary CFT. 
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Figure 6: Energy density divided by iV^/(27r^) as a function of the Schwarzschild radial coordinate 
R. The energy density is finite at the horizon and decays like far from it. 

In Fig. 6 we plot Rq (T^^) for the fixed point of the flow, and in Fig. 7 we show the 
function t4{R) that can be extracted from the numerical solution. We observe that -^(T^"^) 
is finite at the horizon and decays like R~^ far from the black hole. In fact the large R 
behaviour of the stress tensor can be conjectured analytically as follows; in the context 
of the Randall-Sundrum single braneworld model, the linearised gravitational field created 
by a point-like source (of mass M = Rq/2) has been calculated in [41,42]. According to 
AdS/CFT from the knowledge of the gravitational field on the brane we can compute the 
expectation value of the stress tensor of the CFT at strong coupling and large N^. from the 
relation [29]: 

G^, = 167rG4(r^.), (4.19) 
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where we note that we have included an extra factor of two on the righthand side than 
appeared in [29] which using our conventions we expect to be present. Then substituting 
the Hnear calculation of the perturbation to G^y for the point source yields 

which we observe is precisely the asymptotic behaviour of our stress tensor. 

0.05 [p ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ^ 

0.00 ; : 

-0.05 - - 
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Figure 7: t^^R) calculated from our numerical solution by fitting the near boundary region to the 
behaviour in the large R region, t4^{R) decays as — _Ro/(4-R). 

Even though the functions t4{R) and S4{R) are finite on the event horizon this is not 
sufficient to guarantee that (4.18) is regular there. However, since our bulk classical solution 
is a smooth geometry in its interior we expect any geometric data extracted will also be 
smooth, although we note that this is not guaranteed (see for example [43]). We may 
explicitly check the regularity of the stress tensor by returning to Lorentzian signature 
and transforming (4.18) to coordinates which are regular on the (future or past) event 
horizon, e.g., ingoing or outgoing Eddington-Finkelstein coordinates. Then the necessary 
and sufficient condition for regularity is that (T^^) = {T^^) at R = Rq, but in our case 
this can be seen to be consequence of the conservation of the stress tensor, which in turn 
is implied by the bulk equations of motion. Therefore, we conclude that (4.18) is indeed 
regular on both the future and the past event horizon. 

5. Summary and discussion 

In the first part of this paper we have shown that the existence of solutions to the Ricci 
soliton equation (2.1) (with a non-positive cosmological constant) on an m-dimensional 
Riemannian manifold with a boundary is governed by a simple maximum principle. In 
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particular, if the vector field ^ vanishes on all boundaries, this rules out the existence of 
non-trivial Ricci solitons. Following [1] we have considered static solutions to the Einstein- 
DeTurck equation by Euclidean continuation to a Riemannian problem which is elliptic. We 
have extended previous discussion by showing how to implement boundary conditions for 
the Einstein-DeTurck equation in a variety of cases: boundaries with 'modified Dirichlet' 
data (Anderson's fixed trace of extrinsic curvature and conformal class for the induced 
metric [8]), 'mixed Neumann-Dirichlet' data (extrinsic curvature proportional to induced 
metric), and ends where the analytically continued Lorentzian metric is asymptotically flat, 
Kaluza-Klein, locally AdS or has an extremal horizon. Furthermore, using our maximum 
principle, we have argued that for solutions with these boundary types, the existence of 
non-trivial Ricci solitons is actually riilcd out.^^ Thus if one can solve the Einstein-DeTurck 
equation, for example using Ricci-DcTurck flow, one is guaranteed that the solution will 
solve the vacuum Einstein equations (with a non-positive cosmological constant). 

We also considered whether Ricci-DeTurck flow preserves our various boundary and 
asymptotic conditions. We have proved that the flow does preserve the class of static 
asymptotically extremal manifolds - we defined such manifolds as ones whose associated 
static Lorentzian spacctimc contains a smooth static extremal Killing horizon whose near- 
horizon geometry solves the Einstein equation. 

In the second part of the paper we focussed on a particularly elegant application of 
these more formal results. We have numerically constructed a 5d Einstein metric with 
negative cosmological term, whose static Lorentzian continuation is asymptotically AdS 
with boundary metric conformal to the 4d Schwarzschild black hole, and whose metric far 
from this boundary tends to an extremal horizon whose near horizon geometry is that of 
the AdS Poincare horizon. We expect this Einstein metric gives the dominant classical 
bulk saddle point solution for AdS-CFT where the dual strongly coupled CFT is put on 
a Schwarzschild background in the Unruh or Boulware vacua, so that asymptotically the 
theory is in vacuum. Our analysis of the maximum principle and boundary conditions in 
the first part of the paper shows not only how to formulate the elliptic Einstein-DeTurck 
equation and boundary conditions, but also gives the beautiful result that the maximum 
principle implies any solution to this DeTurck equation must be Einstein and not a Ricci 
soliton. 

In situations where one tries to find an Einstein solution that does not exist there are 
two obvious ways in which the methods above can fail. Firstly the solution found can 
be a soliton. Secondly, for any initial data Ricci flow develops a finite time singularity 
that may indicate a solution exists but on a manifold of a different topology. For our 
application it is not a priori obvious that a solution with the prescribed boundary conditions 
docs exist, and thus ruling out the possibility of a soliton analytically is a very nice result. 
Of course, one can always check numerically whether one has found a soliton or not, but it 

'^^With the exception of the mixed Neumann-Dirichlet condition for 'negative tension'. 

^'^This would be the case if one tries to find a positively curved Kahler-Einstein metric on the second del 
Pezzo surface, since such a metric does not exist. A Kahler Ricci soliton does though, and this has indeed 
been found by Ricci flow [7] . 

^®This feature can be seen in Ricci flows relating to static black holes in cavities [38]. 
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is obviously preferable to have analytic control over their existence. If one knows solitons 
cannot exist, and so the DeTurck vector ^ should vanish, then the norm of ^ provides a 
good estimator of numerical error for a solution. 

Another very attractive feature of our example application is that it is particularly 
simple to construct the solution. Whilst black hole solutions often have Euclidean negative 
modes, our solution does not, and any putative mode is projected out as the conformal 
boundary metric is not dynamical but is fixed to be Schwarzschild. Hence we may find the 
Einstein solution very straightforwardly by simulating the parabolic Ricci-DeTurck flow 
with some initial guess. With no negative modes, the solution is a stable attractor, and 
in principle we may flow arbitrarily close to it if we flow for sufficient time. Unlike other 
solutions studied with these methods, such as in [1] there is no critical behaviour to be 
tuned in the Ricci-DeTurck flow, and no need to implement the more complicated Newton 
method algorithm. 

We have constructed the solution to very high accuracy using pseudo-spectral meth- 
ods and have exhibited various properties of the solution. In particular we have found 
the shape of the horizon is such that we expect no Gregory-Laflamme instability to exist. 
Together with the lack of Euclidean negative modes^^ (at least in the static axisymmetric 
class as the solution is a stable fixed point of Ricci-DeTurck flow) wc expect that the solu- 
tion is dynamically stable, although we emphasize that a proper Lorentzian perturbation 
calculation would be required to determine this beyond our heuristic discussion. In partic- 
ular it is important to include the possibility of instabilities involving the internal space in 
AdS/CFT (for example an S^) and localized near the 'tip' of the bulk horizon. 

We have claimed that the IR boundary conditions chosen are appropriate to describe 
the dominant classical bulk saddle point dual to the CFT in the Unruh and Boulware vacua. 
We will conclude the paper with some discussion on the interpretation of the computed 
CFT stress tensor. 

An important point is that for any classical gravity calculation of (Tij) for the dual 
CFT, only the leading 0{N^) behaviour can be extracted, whilst 1-loop bulk gravity correc- 
tions would be required to compute the remaining 0(1) contributions to the stress tensor. 
Thus we emphasize that wc have not computed the full stress tensor for the CFT but only 
the leading O(N^) part which is common to both the Unruh or Boulware vacua, these 
being distinguished only at 0(1). 

Already at this leading order one can distinguish the Hartle-Hawking vacuum from 
the Unruh or Boulware ones, as for these latter cases we expect ■^{T,-') ^ as i? ^ oo, 
whereas for the Hartle-Hawking vacuum one expects the energy density is 0{N^) as R ^ 
oo. In the gravity solution, in order to describe the dominant classical bulk saddle point 
for the CFT in the Hartle-Hawking vacuum, one would take a finite temperature horizon 
in the IR, at the same temperature as the UV horizon, rather than the extremal horizon of 

^■^Thcre is of course no direct relationship between a Euclidean negative mode and a Lorentzian unstable 
mode as one cannot simply continue a time dependent mode and its boundary conditions between the 
two signatures, so a correspondence only holds for static perturbations [44]. That said, the existence of 
a Euclidean negative mode (other than that associated to local thermodynamical instability ic. a Gross- 
Perry- Yaffe mode, which does not occur here) might be taken as a weak indication of a potential instability. 
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Poincare AdS [13]. We note that our solution docs describe some classical bulk saddle point 
dual to the CFT in the Hartlc-Hawking vacuum, but it is not expected to be the dominant 
one.^° From a Euclidean perspective this can be understood by considering which classical 
saddle dominates the free energy, but also from the expectation of strong corrections to 
the semiclassical saddle point in the extremal region due to degeneration of the Euclidean 
time circle which is periodically identified (with antiperiodic fermion boundary conditions) 
due to the finite temperature boundary black hole horizon [45]. 

The leading 0{N^) behaviour of our stress tensor is very interesting as since it is 
static, there is no energy flux to infinity and the stress tensor is regular on the past and 
future horizons. Our expectation is that by including 1-loop graviton corrections in the 
bulk which will determine the 0(1) part of the stress tensor at large Nc, the stress tensor 
will be regular only at the future horizon in the Unruh case, and will be singular on 
both past and future horizons in the Boulware case as is usually found. What makes this 
picture so interesting is that if the CFT were a free theory then on rather general grounds 
(for example [46]) taking the Unruh or Boulware vacua already would yield the leading 
O(iV^) part of the stress tensor to be singular on the past horizon (Unruh), or past and 
future horizons (Boulware), and in the Unruh case one would also sec a flux of radiation 
in the stress tensor out to inflnity. Thus the strong interactions of the CFT appear to 
significantly change the behaviour of the theory from that expected by intuition from free 
field theory. We note that this is very much related to the issue of existence of Randall- 
Sundrum braneworld static black holes [13,28-33], where we believe the arguments against 
their existence based on extrapolation from free field theory would also incorrectly imply 
non-existence of the solution we have found here. 

We now present a simple physical picture which we believe captures the physics of 
this strongly coupled CFT in the Unruh vacuum. The strong self-interaction of the CFT 
is attractive (corresponding to the attractive 5d gravity force in the gravity dual). Thus 
the thermal radiation pressure from the black hole horizon must compete with the strong 
attractive self-interaction of the resulting plasma, which naively wishes to collapse back into 
the horizon of the fixed background, thus creating a "halo" of plasma near the horizon. 
It appears that at order 0{N^) in the stress tensor a static equilibrium is reached, with 
a thermal halo of plasma forming, and at every point in the halo the radiation pressure 
balances this long range attractive self-interaction.^^ Although the halo has no edge and 
extends all the way out to infinity, since -^{T^'^) = 0{R~^) at large R, the energy of the 

c 

halo is finite. While at 0{N^) static equilibrium is reached and there is no energy flux, as 
we have mentioned above, the 0(1) component of the stress tensor will presumably exhibit 
Hawking radiation out to infinity in the Unruh vacuum. 

It is an interesting question whether the static 0{N^) halo will persist to finite 't Hooft 

^°This is analogous to considering the CFT in flat space at flnite temperature where the Poincare-AdS 
solution does give a bulk classical saddle point but is dominated in the bulk partition function by the planar 
AdS-Schwarzschild saddle point. 

^^Note that there is also the gravitational attraction to the horizon due to the background metric, but 
this is a small effect, and would not be expected to be strong enough to confine a thermal halo for a free 
CFT. 
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coupling, or is only an artifact of being at infinite 't Hooft coupling. In principle one might 
imagine studying this by including the leading a! corrections in the gravity dual. 

Another interesting direction is to study the case of the Hartle-Hawking vacuum, which 
has been considered in [13,14] where solutions are conjectured to exist and have been termed 
black 'droplets' and 'funnels'. In the case of droplets, the construction of solutions would 
follow essentially what we have done here, the only difference being that the extremal 
horizon boundary would be replaced by a non-extremal horizon 'fictitious' boundary. Our 
maximum principle arguments would again rule out the existence of solitons, and further- 
more solutions would be expected to be stable fixed points of Ricci-DeTurck flow. 
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A. Static extremal Killing horizons 
A.l Coordinate systems 

In this section we prove that any static spacetime containing a smooth static extremal 
Killing horizon, can be written in coordinates (t, x") valid outside the horizon > 0, 
such that the metric takes the form (2.28) where (2.29) and (2.30) are satisfied. Our only 
assumptions are that the static Killing vector field is timelike just outside the horizon and 
that cross-sections of the horizon are simply connected. 

We begin by noting that the metric in a neighbourhood of any extremal Killing horizon 
of a Killing field V can be written in Gaussian null coordinates (f , r, x") as: 

g = 2dv ^dr + r ha{r, x)dx"' — -r'^F{r, x)dv^ + ^abi^, x)dx°'dx^ (A.l) 

where V = d/dv, the horizon is r = 0, and are coordinates on cross-sections of the 
horizon which we denote by Ti. Note that in these coordinates all metric functions are 
smooth at r = 0. 

We are interested in static extremal horizons, i.e. when V is hypersurface orthogonal, 
so y A dV = everywhere. It can be shown that V is hypersurface orthogonal iff 

dh = rhA drh (A. 2) 

dF -hF + r{Fdrh - hdrF) = (A.3) 

where d denotes the exterior derivative at constant r, so e.g. dh = d^ahbjdx"" A dx^, and 
drh = (drha) dx'^ . By smoothness we may write 

F = Fo(x)+rFi(x) + 0(r2) , h = ho{x) + rhi{x) + 0{r^) (A.4) 



-38- 



where FQ,Fi,hQ,hi are smooth on Ti. Substituting these into (A. 2) and (A. 3) gives 



dho = 0, dFo = Foho , (A.5) 

dhi = hoAhi , dFi = 2Fi/io . (A.6) 

The two equations (A.5) are easily solved 

ho = dX{x) , Fo = coe^(^) (A.7) 

where cq is a constant and A(x) a function on T-L, which agrees with the general form for 
a static near- horizon geometry found in [11]. For simply connected Ti these expressions 
are valid globally. The assumption that V is timelike just outside the horizon, implies 
Co > 0, and w.l.o.g. we will set cq = 1. The near-horizon geometry can then be written as 
(2.27) where g = e'^r and Tq = e~^. We are now interested in generalising this to the full 
spacetime metric. For later convenience we note that we can solve (A.6) 

= dXi{x) , Fi = cie2^(^) (A.8) 

where Ai(x) is a function on H and ci a constant. 

It is clear one could proceed by solving (A. 2) and (A. 3) order by order in r. Instead 
we now find new coordinates which make the static isometry manifest. Assuming r > 
and F > (as is the case if V is timelike just outside the horizon), it can be shown that 
V A dV = implies there is a coordinate transformation defined by 



t = v-f{r,x), df = ^— (A.9) 



such that the spacetime metric (A.l) becomes manifestly static 

g = -r'^Fdt^ + r'^Fdf + -iabdx^dx^ 
] / dr \ 

= -r'^Fdf + - — + hadx" + -fabdx'^dx'' . (A.IO) 
i* \ r J 

It is convenient to change radial coordinate aswell. Thus define 

g = Tr (A.ll) 

for some positive smooth function F chosen such that F|g=o = e^^^\ The metric then 
becomes ^ 

g = -g^Tdt^ + r(^^ + QUJadx"^ + -fabdx'^dx'' (A.12) 

where 

T=Y2, R=-\l ^1 , ha= — + gn f^)'^'^' 

define the smooth functions T,R and 1-form oja- Note that as defined in equation (A.13), 
LOa is smooth at = as a consequence of the fact that Hq = d\ and the choice Fq = e'^^^^. 
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We have thus derived the form of the metric claimed in the main text (2.28). Fur- 
thermore, since we have related this to Gaussian null coordinates we may now deduce 
conditions on T, R, oja, "fab which are necessary and sufficient for the existence of a smooth 
static extremal horizon. First it is clear that T\g=Q = R\g = which confirms (2.29). 
Next, defining Ti = dQT\g=Q and Ri = dgR\Q=Q, one can show that 

Ti-i?i = ^^ (A.14) 

where Fq,Fi are defined by (A. 4). Now, using the equations for Fq,Fi worked out above 
(A. 7) and (A.8) implies 

^ = 2c, (A.15) 

which confirms (2.30). In the main text it was shown that given the metric (A. 12), these 
conditions are sufficient to guarantee a smooth extremal horizon at g = 0. 

This completes the proof that outside any smooth static extremal Killing horizon, the 
metric can be written as (2.28), where T, R, uJa, Jab are smooth functions at = together 
with (2.29) and (2.30), assuming the static Killing field is timelike just outside the horizon 
and Ti is simply connected. 

A. 2 Ricci-DeTurck flow of asymptotically extremal manifolds 

In this section we show that Ricci-DeTurck fiow preserves static asymptotically extremal 
Riemannian manifolds. These were defined in the main text and are equivalent to the 
class of Lorentzian spacetimes containing a smooth extremal static Killing horizon. As 
argued, we may choose coordinates such that near the horizon this class of metrics can be 
written as (2.34). We already argued that the near-horizon geometry is preserved, simply 
as a consequence of the fact it is a fixed point of the fiow equation. Now, we derive flow 
equations for the "first-order" quantities Ti = dgT\g=Q and Ri = dgR\g=o. The obvious 
thing to do is calculate the Ricci tensor for the full metric (2.34) and Taylor expanding to 
first order. However, this is somewhat cumbersome, and we will take a different approach. 
Instead we will treat the first order terms as a linearised perturbation about the near- 
horizon geometry. 

Set g = er and t = t/e for e > 0. The full metric (2.34) is then: 

,^T..^.*^.fl,(^ + ™y.")%,.,..W ,A.10, 

where = T{er,x) etc. This is a 1-parameter family of metrics which coincide with gj\ff{ 
when e — >• 0, i.e. ge ^ 9o = 9nh as e — )■ 0. Now define the first order perturbation by 



, dge 



(A.17) 



Similarly one can define the 1-parameter family of reference metrics g^ , which also satisfies 
ffo = 9nHi and h = {dg / de)e=o- Differentiating the full fiow equation dxg^i, = — 2(i?^;^ — 
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^(fi^i/) ~ ^9fiu) wrt e implies a flow equation for h^i, 

dh 



22 



(A.18) 



where = V^h^^i — \d^h, A^, is the Lichnerowicz operator, and all quantities are with 
respect to the "background" metric go = g^H (2.27). Explicitly, the Lichnerowicz operator 
is given by 

Al V = -V V + 2^P(M V - '^Ri^pu.hf" (A.19) 



where V and and R^upa are the Levi-Civita connection, the Ricci tensor and Riemann 
tensor of gjyH- Therefore all curvature calculations are with respect to gjyn (2-27), making 
the task considerably simpler. In our case one can check 



h = Tir^dt^ + 



+ 2Toujldx°-dr + r-f^t,dx°-dx 



1 J^a^h 



(A.20) 



where Ti = dgT\g=o, Ri = dgR\g=o, uj^ = u}a\(,=o and 7^^ = d(,-/ab\e=o all depend only on 
X. Since we are interested in obtaining an evolution equation for Ti — Ri we need only 
consider the tt and rr components of (A.18). 

Working in the coordinate basis defined by x^ = (t, r, x") gives the following non-zero 
Christoffel symbols for the metric g^n (2.27): 



' ^tt — ' ^ rr 



-2-pa _ „2-pa 
' ^tt — ' ^ rr 



1 



nc = r[7' 



1 



■ tr 



ta 



dgTo 
2Ta 



Oia 
be 



The relevant parts of the curvature tensors are then 



R 



trtr 



r-^Rt. 



-To 



r ^Rtatb = r^Rrarb 



-VbdaTo 



dgTodbTp 
2To 



1 



V^Tn - 1 



(A.21) 



,(A.22) 
(A.23) 



where V and • now refer to the Levi-Civita connection and contraction with respect to the 
metric 7^^^, respectively. It can then be shown that 



Ai/itt 



Ar/l, 



L'^rr 



2 dTo ■ dTi _2J\_ {Ti-Ri){dTof 
To To 2Ti 

daTodbTo \ 
To J 

2 dTo • dRi 2Ri (Ti - i?i)(dTo)2 
-V /ii H — — h 



-2dTo ■ Lo^ + 7f ( VaQfeTc 



To 



To 



2T^ 



-2dTo ■ LO^ + 7f ( Va^feTc 



daTodbTo\ 
To J 



(A.24) 



(A.25) 



^^We have used R^uA^=o = f AL/ip,y + V(^u„) where = V/i^^ - and ■£ (2V(,iCi/)) = ^io^^"' + 

2V{^^^)|e=o, and = ?e=o = and ^^1^=0 = " ^m- Note we are using go = ^o- 
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Therefore we deduce 

-^ALhtt - rAihrr = -V\Ti - + — • d(Ti - Ri) . 

(A.26) 

We also need VtVt and VrVr- It is easily checked that vt = and drVr = 0. This is sufficient 
to show that 



-^1 J- \/ 71 — 0) J- 

2r2 



tVt = r'^Vr H , VrVr = -Vr H (A.27) 

and therefore 



\vtvt - rVrVr = . (A.28) 



Putting all this together, we have the flow equation 



- ^i) = v^m - - • d{T, - R,) + m - R,) ( ™! + ^ 



Now, define the function 

^ = . (A.30) 

Note the tt component of the near-horizon geometry Einstein equation reduces to 

V^To + 2 + 2ATo = . (A.31) 

Combining this with the flow equation for Ti — Ri (A. 29) one can derive a remarkably 
simple flow equation for ip: 

|^ = VV. + ||^.*. (A.32) 

This shows that if ip is constant at A = 0, then by uniqueness of the flow equation it 
must be constant for A > 0. This is sufficient to establish that asymptotically extremal 
manifolds are preserved by Ricci-DeTurck flow, provided the near-horizon geometry solves 
the Einstein equations. 

B. Static non-extremal Killing horizons 

Consider any static spacetime containing a smooth, static, non-extremal Killing horizon of 
V, such that V is timelike just outside the horizon. In this section we prove that outside 
the horizon the metric can always be written in coordinates {t,w,x'^) such that the metric 
takes the form (2.37) for u; > where tt; = is the horizon. 

We begin by noting that the metric in a neighbourhood of any Killing horizon of a 
Killing field V can be written in Gaussian null coordinates (u,r, x") as: 

g = 2dv ^dr + r ha{r,x)dx"' — -rf{r,x)dv^ + '~^ab{r,x)dx°'dx^ (B-1) 

where V = d/dv, the horizon is r = 0, and x"" are coordinates on cross-sections of the 
horizon which we denote by Ti. Note that in these coordinates all metric functions are 
smooth at r = 0. 
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We are interested in static horizons, i.e. when V is hypersurface orthogonal, so VAdV = 
everywhere. By smoothness we can write / = fo{x) + 0(r) for some smooth function 
/o on Ti (and similarly for ha and Jab)- Expanding the trx" component of the staticity 
condition V A dV = around r = 0, then allows one to show that dafo = 0. Therefore 

/o = 2k (B.2) 

where «; is a constant which is easily shown to be the surface gravity. If k = then the 
horizon is extremal, and therefore in this section we assume k ^ 0. Since we assume V is 
timelike just outside the horizon we deduce n > 0. 

Staticity implies that outside the horizon r > 0, there exists a coordinate transforma- 
tion defined by 

dr + rh ^ 

t = v-X , dX = — (B.3 

rf 

in terms of which the metric becomes 

g = -rfdt^ + —Adr + rhadx°-f + -fabdx'^dx'' . (B.4) 
rf 

Now, for r > 0, change the radial variables by defining 

w"^ = r(r, x)r (B.5) 
for some positive smooth function T. Then, for w > 0, the metric becomes 

g = -w'^T^dt'^ + W^{dw + wQadx^'f + -iabdx'^dx^ (B.6) 
where we have defined 

where T,W > 0, which is indeed of the form claimed in the main text (2.37). Notice that 
we have used the assumption that / > near r > 0, i.e. V is timelike just outside the 
horizon, to define T, W. 

We may now deduce necessary and sufficient conditions on the metric functions T, W, Oa, 'jab 
for (B.6) to describe a smooth non-extremal horizon at w = 0. First note that in the orig- 
inal coordinate r these metric functions are smooth at r = 0, and hence we deduce they 
must be smooth in (not w) at w = 0. Also we have 



f 



10=0 



(B.^ 



r=0 



where in the last equality we used (B.2). As we show in §2.2.5, these conditions are sufficient 
to ensure that the metric (B.6) describes a smooth non-extremal horizon at u; = 0. This 
completes the proof of the claim in §2.2.5. 
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C. A comment on Ricci flow in ALH manifolds 



The idea of the Ricci-DeTurck flow, and indeed the Einstein-DeTurck equation is that 
one can avoid choosing an expHcit coordinate system in order to gauge fix the Einstein 
equations. Rather the DeTurck vector provides extra degrees of freedom to ensure the 
equations are parabohc and eUiptic respectively. However, for Ricci flow in ALH manifolds, 
an obvious question is whether the Fefferman-Graham (FG) expansion is valid. Although 
somewhat outside the scope of this paper, for completeness we present a simple argument 
which indicates that this is not the case. 
Consider the scalar diffusion equation 

dx(t> = VV (C.l) 

in the n-dimensional hyperbolic background 

(C.2) 



9 



dp^ hij(x)dx'^dx^ 
4p2 p 



which satisfies Ric{g) = — (n — l)i iff Ric(/i) = 0. In these coordinates the defining 
function = p. It is easy to show that explicitly it is given by 

^ = ^[-2(n - 3)<^' + V2<^] + ^ (C.3) 

where (p = dx^' and (j)' = dpcf) and V is the connection of h. 

For simplicity we will consider only the case of even dimension n. Then we may try to 
find a Frobenius expansion for the scalar flow as one does for the case of Laplace's equation. 
Such an expansion will have the form, 

<t^{P,x) = E (<^2fc(x)/ + 02fc+n-l(x)p("-l)/2+'^) (C.4) 
A:>0 

and we note in odd dimensions one can perform a similar expansion but require logarithmic 
terms. Expanding and solving (C.3) order by order in p gives the following relations: 

00 = , 0„_i = , 

; m, ^ 1^2, r I m > , m even (C5) 

0m = -^oi'^ - 1 - + 0^-2 , for < 

£^ i'^ I m > n — 1 , m odd 

We see that the two terms in the expansion (C.4) are not coupled together by the recursion 
relation. Consider first the integer powers of p. The boundary value of (j) (the 'non- 
normalization data' in AdS-CFT language) given by 0o, does not evolve and we may solve 
for 02 as a function of A: 

ho 



V^0n _2(n-3)A 



where U2 is independent of A. For 2 A; < n — 1 we can recursively solve for 02fc and find that 

V^02(fc-l) 



^2k 



2k{n -l-2k) 



+ ... (C.7) 
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where . . . represent terms which decay exponentially in A, and therefore (p2k converges to 
a value determined by (pQ. Then (j)n-2 is the last coefficient which decays exponentially in 
A, and for 2k > n we see that (p2k a^^e still determined however they grow exponentially 
in A and hence do not converge at late flow time. For the half integer powers the leading 
term 4>n-i (the 'normalization data') does not evolve. However all the higher odd terms 
4'2k+i grow exponentially in A, and so again do not converge. 

This example shows that the small p expansion (C.4) for A > breaks down at 
, since for large m one gets terms in the p expansion which grow exponentially 
as ~ exp(m2A/£2)/9™'/2. One obtains a similar result for odd dimensions. 

We expect analogous results for the Ricci flow thus indicating that the FG expansion 
in p cannot be used to describe the flow in some constant neighbourhood of the boundary. 
This means that there it is unclear how to define the boundary stress tensor along the 
flow. We note though that, curiously, the lower order terms 2A; < n — 1 in the expansion 
all converge exponentially to the values at a fixed point (where the FG expansion is valid), 
just as in the example of the scalar field above. 
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